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ABSTRACT 


In the past decade there has been renewed scientific interest 
in Extremely-Low-Frequency (ELF) radio waves, in the range from 3 Hz 
to 3 kHz. A brief historical review is included of investigations 
related to ELF phenomena, and the spectrum characteristics are out- 
lined with special reference to the Earth-ionosphere cavity "Schumann 
Resonances" and to transient disturbances related to global lightning 
activity. The current technology of measurement at ELF frequencies is 
briefly examined and the design criteria are discussed for a special 
field-portable receiver and magnetic-tape recording system for ELF 
sienals. The detailed design equations and schematics of a tested 
prototype system are then presented, with special emphasis on the sub- 
systems involved: a convenient charge-induction type of E-field sensor, 
a modified 60-Hz notch filter with a controlled transient response, a 
Pulse-Duration—Modulation (PDM) method for recording the ELF signals 
on magnetic tape, and a PDM demodulator system with a tape-speed- 
variation compensator, to recover the ELF signals from the magnetic 
tape. A number of examples are included of field-recorded natural ELF 


transients and background noise. 
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CHAPTER 1 


INTRODUCTION 


1.1 THE EXTREMELY-LOW-FREQUENCY (ELF) SPECTRUM 
1.1.1 HISTORICAL BACKGROUND 


The electromagnetic spectrum above 100 kHz has been extensively 
exploited for many years. The commercial incentives have ensured adeq- 
uate theoretical and experimental investigations of propagation phenom- 
ena and the continuing development of communications systems. 

Below 100kHz, the greatly increased size and cost of transmitt- 
ing antennae, the decreasing radiation efficiency and the limited chann- 
el capacity have almost eliminated commercial interest. However, in the 
Very-Low-Frequency (VLF) range, nominally 3 kHz~-30kHz, very reliable long- 
range navigation and time-standard systems have been developed [1],[2], 
eile herattenuatton raterrs) less than 10° ¢d8/1000km at VLY5 and 
propagation is primarily via waveguide modes in the earth-ionosphere 
cavity, so that such systems are relatively immune to the normal vagar- 
ies related to ionospheric disturbances. 

Below these frequencies again, at Extremely—Low—Frequencies 
(ELF), nominally 3 Hz-3 kHz, lies a portion of the electromagnetic spectrum 
which in recent years is being intensively re-examined. 

Historically, about 1900 Nikola Tesla observed periodic excit-— 
ation of tuned spark-gaps by receding lightning storms, and envisioned 
broadcast power based on the resonant pumping of the earth's charge at 


about 10Hz [5]. His proposals for world-wide communications on the same 
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theme [6] were never developed, and after Marconi's successful trans- 
atlantic signal experiments in 1901, commercial exploitation concent- 
rated on the higher frequencies. Barkhausen reported in 1919 [7] that 
sharp audio whistles were often noted in the high-gain amplifiers used 
to monitor enemy telephone circuits (via earth-return currents) and in 
1930 he postulated that these whistles might be due to a long-path 
ionospheric dispersion of the audio-frequency components of far-distant 
lightning flashes [8]. 

Some theoretical justification for such a dispersion mechanism 
was offered in 1935 by Eckersley [9] but it was not until 1954 that 
Storey [10] offered the now accepted theory that ionospheric disturb- 
ances due to lightning flashes can find propagation paths far beyond the 
ionosphere along the Earth's magnetic-flux lines, so that inter-hemis- 
pheric transmissions and reflections are possible. Extensive investigat- 
ion of these 'whistlers' has been made since that time because of their 
relation to the characteristics of the ionosphere and the magnetosphere 
Pi) s 

The possibility of earth-ionosphere cavity resonances again 
became of interest in 1952, when Schumann offered a mathematical analysis 
ofthe subject [12])[13] and showed that for the ideal case of perfectly 
conducting earth-ionosphere boundaries, the eigenfrequencies would in 
fact Heait \.59/..46 A Ge UL za en". one OG zi, 08 23:02, 25:00) ah ee) . 
The general background noise at ELF makes it difficult to directly obse- 
rve one transiently excited resonances, although Schumann and Konig 
in 1954 [14] offered several recordings of damped oscillations at about 


9Hz showing field strengths in the order of 500 microvolts/metre. 
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It was not until 1960 that Balser and Wagner [15] secured 
averaged power spectra of ELF noise levels which clearly verified the 
existence of the Schumann resonances. The measured eigenfrequencies 
Were auuiaudizenlasuizee2O7oeZ 12) Shand were sufficiently xesolved 
so that estimates of the cavity Q could be made from the lowered 
frequencies of the mode peaks and their relative amplitudes. Subsequent 
verifications of these spectra have been made by other investigators 
using different methods, and military interest has led to much experim- 
ental and theoretical work on ELF propagation characteristics and 
ionospheric conductivity profiles. Extensive bibliographies of relevant 
ELF investigations have been included in papers by Galejs [16], Jones 
iA) waWaut Bhs | Bandtothens: 

In the past decade particularly, the major research into ELF 
phenomena has been funded by the U.S.Navy under ‘Project Sanguine". 
Even one-way communication with submerged submarines on a global scale 
is considered imperative and appears feasible; below 100Hz attenuation 
is less than 1 dB/1000Km and sea-water penetration can be in the order 
of tens-of-metres. Unfortunately, transmitting antennae must cover 
literally thousands of square miles and require megawatts of input power 
for a few watts radiated. For one case cited, the radiation resistance 
for an ELF antenna 620 miles lone was only about 0-Ol ohms, so that 
less than 0.03% of the input power was radiated [19]. (These studies 
have been of such wide interest that a special issue of the IEEE Trans— 
cone on Communications, Vol. COM-22, No.4, April 1974, was devoted 


entirely to papers related to Project Sanguine.) 
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Because of the very high ground currents associated with such 
an ELF transmitting array, and their distribution over such an immense 
area, the potentially adverse environmental and biological effects 
were extensively investigated [20] but there remains a wide-spread 


controversy over system cost-effectiveness and environmental impact [21]. 


tii. 2 SPECTRUM CHARACTERISTICS 


Since the ELF wavelengths are much greater than the height of 
the ionosphere, propagation between the earth-ionosphere sheet conduct- 
ors is essentially in the TEM mode. However, the boundaries do have 
finite conductivities resulting in small horizontal E-field components, 
and these lead to wave leakages into the bounding media. In practice, 
this energy loss is a very small fraction of the propagating wave energy, 
so that the attenuation rate is very low. 

Propagation characteristics at ELF over various path lengths and 
under a wide range of conditions have been measured by many investigators 
[22],[23],[24], giving the typical attenuation curves shown below . (The 
dip in the region of 10kHz is due to a resonance absorption in the earth- 


ionosphere cavity.) 
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Figure 1.1 ELF Attenuation vs Frequency (from Maxwell [33]) 
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The natural ELF noise spectrum consists primarily of random 
disturbances related to world-wide lightning activity. The wide spectrum 
of electromagnetic radiation generated by lightning strokes has been 
earerully snvestieated for many years [25],[26],[27],1(28],[29] with 
various portions of the radiation spectrum attributed to specific flash 
processes. Most of the energy radiated in the 30 kHz - 100 kHz range is 
related to the very fast pre-discharge "leader" breakdown processes, 
and attenuation with distance is sufficient to limit the propagation 
to a few thousand kilometers. The main-discharge "return-stroke" 
radiation peaks at about 10 kHz, while slower decay currents between 
the multiple return strokes (described as "slow tails") radiate 
considerable energy below 1 kHz. At such extremely-low frequencies, the 
very-low attenuation rate permits world-wide propagation, and this 
COnsteiliures tenes hbackeround noise., (30), 13h) 5132151331. 

The ELF noise below 1 kHz may actually be considered as being 
made up of two components: 

1) A background level of impulse noise, with a Poisson repetition- 
rate distribution centred on a world-wide lightning-flash 
average of perhaps 100/sec, with a typical power spectral 
density of about 0.05 (mV/m)2/Hz [32]. 

2) A superimposed random noise, with peak-signal levels one or 
two orders of magnitude greater than the background (but with 
a much lower average repetition rate) generated by local or 
regional lightning storms. 

Ties resultant ELF noise is typically as shown in Figure 1.2, which 


shows narrow-band (5 Hz-40 Hz) recordings of the horizontal magnetic-field 
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components, recorded simultaneously at widely-separated geographically- 
located points. The time correlations of many of the noise bursts are 


quite evident. 
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Fig. 1.2 Narrow Band ELF Noise Records (from Shand [34]) 


(Horizontal Magnetic-Field Components , 5 Hz-40 Hz) 


These noise bursts are superpositions of resonant-mode disturbances, 
with only the first two modes being evident in these recordings. Diurnal 
variations in activity produce changes in the number of bursts which 
occur, but in general the average amplitudes remain essentially the 


same. 
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Wide-band receivers can of course display the impulsive 


nature of ELF signals with less modification by the receiving system: 
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ELF wide-band (20-300-Hz) noise. (a) Saipan median level. (b) Saipan high level. (c) Norway. 


Fig. 1.3 Wide-Band ELF Noise (from Bernstein [35]) 


(recording methods and scaling not specified) 
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Occasionally, high-speed chart recordings may display the 


sequence of pulses generated by the multiple strokes of a lightning 


flash: 
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Examples of a pulse train of ‘slow tail’ ELI atmospherics recorded simultaneously at 


the Jlawan and Arizona Stations. 


Figure 1.4 Lightning-Flash Pulse Trains (from Hughes [36]) 


Despite the continuous presence of such ELF activity, well- 


defined oscillations at Schumann-mode frequencies are difficult to obse- 


rve. Typically, oscillatory disturbances retain their phase coherency 


for only a few cycles before the random superposition of another disturb- 


ance. An isolated mode oscillation such as that shown in Figure 1.5 is 


therefore a comparatively rare occurrence, 


Nevertheless, Shand [34] 


suggests that "excepting the sources of local lightning and man-made 
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interference, no measurable non-Schumann activity has been found in the 


frequency range 5 - 40 c.p.s.". 
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Figure 1.5 A Schumann-Resonance Fundamental-Mode Transient 


(fromeCeawanetcdlen 15710) 


Due to the difficulty in measuring the Schumann resonances 
directly, most of the ionospheric research based on the changing 
characteristics of the resonances has relied on averages of power 
spectra integrated over many minutes. (Short-period spectral estimates, 
as shown in Figure 1.6, may show considerable fine structure and even 
Peak splitting, but, Rycrorte [32] considers Such detail to be spurious, 
due to statistical fluctuations in signal strensth and finite record 


lengths.) 
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Figure 1.6 
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Short-period Spectrum Measurements 
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Short-Period ELF Power Spectra 
(from Larsen and Egeland [38]) 
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2 EXISTING ELF MEASUREMENT TECHNIQUES 
Mena H-FIELD SENSORS 


At ELF, induction coils are of course very inefficient due 
to the extremely low rate-of-change-of-flux of the magnetic-field 
components. For useful levels of induced signal voltages, such inductors 
therefore tend to be a problem because of their physical size, and from a 
stability standpoint. 

The system used by the Pacific Naval Laboratory (Victoria, B.C.) 
is typical [39]. For the horizontal H-field components, electrostatically- 
shielded inductors are used, each with some 30,000 turns of #22 wire 
wound on a laminated mu-metal core some 72 inches long, with a cross- 
section of approximately 1.4 sq.in., and a nominal weight of about 
94 pounds. Inductance is given as 600 Henries, and sensitivity is reported 
to be approximately 150 microvolts/gamma/Hz. These coils are rigidly 
supported and shielded from wind buffeting, since mechanical vibration 
of the coil in the earth's magnetic field can easily generate signal 
artifacts an order of magnitude (or more) greater than the ELF signals 
of interest. 

For the still weaker traces of vertical H-field components, 
Barkhausen noise in the high-permeability core becomes unacceptable, so 
that very large air-core loops must be used ( typically a buried shielded 
Loop,oreseto, LO turns, ,with a,doopediameter,of L000. ft. or more Ne 
Associated receivers and recording equipment are typically at least 
500 ft. away from the sensing coils, with connections made via buried 


shielded cables. 
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For geophysical applications, where the ELF signals are integ- 
rated over 30 seconds or more, with a number of such readings averaged, 
considerably smaller coils can be used, but vibration in the Earth's 


magnetic field remains a problem [40],[41]. 


Loe E-FILELD SENSORS 


As noted previously, far-field propagation at ELF is via a 
TEM mode, so that above ground, only the vertical E-field component is 
normally measurable. Also, since the height of a practical E-field 
antenna will be a very small fraction of a wavelength at ELF, its 
effective source impedance will be very high and almost totally capac- 
itive. 

The original work on the Schumann resonances by Balser and 
Wagner at M.1.7..[45] made use of a 120-ft. metal tower as the vertical 
E-field sensor,’ while Rycroft at Cambridge, England [32] used a 190- 
metre horizontal wire suspended at an average height of 4.5 metres 
above the ground. Much shorter whip antennae have been used for special 
applications; Maxwell made use of an 8.85 metre vertical-whip antenna 
for ELF noise-field measurements [42] and Hughes used a 21-ft. vertical 
antenna for making field recordings of ELF lightning-flash waveforms 
L437 

The absolute calibration of such vertical-whip antennae can be 
rather difficult, since the “effective height" of the system is a funct- 
ion of antenna capacitance to ground, all mounting and amplifier input 
capacitances, ground-plane effectiveness, and system geometry. The 


method suggested by Clayton et al. [44] made use of a near-by reference 
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monopole to generate a calculated [-field, and a special Schering 
capacitance bridge to measure the effective whip-antenna capacitance. 

A practical problem related to the vertical-whip antenna is 
the wind—induced swaying and vibration, which can generate spurious 
signals by modulating antenna capacitance. The low-frequency artifacts 
below 5 Hzcan be minimized by a controlled low-frequency roll-off in 
the following amplifiers, but any vibration components above 5 Hzare 
directly in the frequency band of interest. 

Two approaches have been used to avoid such motion artifacts; 
a tripod made up of three 6-ft. aluminum tubes ( with the preamplifier 
in one leg ) has been used by Jones and Kemp [45], and a ball antenna 
rigidly mounted on a supporting mast some 15 metres above a building 
roof has been used with success by Ogawa et al. in Japan [46]. 

Horizontal E-field components at ELF appear mainly as leakage 
fields in the conducting boundaries, and are related to the local earth 
conductivity. Such horizontal-field components are usually measured in 
terms of the voltage difference between widely-spaced ground electrodes. 
The voltage measured by such a grounded horizontal-wire antenna is 
related to the above-ground vertical E-field, to a good approximation, 


by the expression [47] 


V=7.45 x 10 °E, 2 fE cosé 
oO 


(ie ee eG Gan E =10mV/m, and antenna length (2%) of 100 m, a frequency 
Ge) Gr UOHZs aesoile conductivity .(o) "Or 10-mhos/n, and o=0 for a wave 


propagating lengthwise along the antenna, V ~ 700 uvolts.) 
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While such signals appear to be of sufficient magnitude to be useful, 
they tend to be obscured by electrode noises and by artifacts traceable 
fomerounu viprarvonlvor long cables. As early as 1945, Dahlbere [48] 
used 'porous-pot" ground electrodes a few hundred feet apart and measured 
signals in the 0.1Hzto10Hz range with amplitudes of a few hundred 
microvolts, but because of the interfering noises they appeared to be 
of little use for geophysical studies. 

More recently, in 1969, Slankis et al. [49],[50] have used the 
8 Hz Schumann fundamental-mode signals for certain types of geophysical 
surveys, using only simple 3/8 in. diameter brass-rod electrodes with 
about 100-ft. spacing. They report integrated average telluric-field 
strengths of approximately 10 my/Km in a sharply-tuned band at 8 Hz. 
(Integrated average-noise amplitudes for periods of 40 sec. to 60 sec. 
were used, and no attempt was made to observe characteristic ELF wave- 
forms.) 

For ELF communications with submarines, a towed horizontal- 
wire antenna in sea water has been used, with the exposed electrodes 
made of titanium, and with electrodes some 1000 ft. apart [51]. The 
high electrical conductivity of the sea water reduces the available 
signal voltages to a very low level, so that special signal encoding 
and extensive digital data-processing of the received signals is 


necessary for a data rate of even l bit/sec. 
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CHAPTER 2 


GENERAL DESIGN CONSIDERATIONS 


From the brief review offered in the preceding chapter, it is 
evident that most of the research into ELF phenomena has required the 
support of rather extensive instrumentation systems. Some of these 
systems required thousands of feet of buried cable and well-equipped 
instrumentation vans, while others were designed around permanent 
elaborate facilities and were part of an international monitoring 
network. Such systems are typically concerned with major and long- 
established research programs or have been developed under military 
research contracts. As a result, independent access to such facilities 
is normally impractical, and even partial duplication can be financ- 
dally prohibitive. 

Nevertheless, many ELF phenomena are of fundamental theoret- 
jealpinterest, so that some ready access to real-time field measure— 
ments is highly desirable for specialized studies as well as for 
general research and teaching purposes. These requirements have led to 
the development of a special ELF recording system, based on the 


following design considerations. 


Dele ePORTABLLITY 


For convenience, the field equipment should be readily port- 
able by a single person. This requirement is related to the fact that 
60-Hz power-line interference lies within the ELF range of interest 


(5Hz - 200Hz) and despite the effectiveness of a 60-Hz notch filter, 
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it is often desirable to make field recordings in selected areas 

as remote as possible from power lines. Since such isolation may also 
be beyond immediate access to motor vehicles, equipment weight is an 
obvious concern. Such portability implies battery operation of amplif- 
iers and recorders, which in any case frees the system from dependence 
on a vehicle-mounted 60-Kz auxiliary power unit, with its attendent 
cost, inconvenience and added electrical interference. 

Some form of magnetic-tape recording is desirable, since this 
permits Eield recording in realtime, with the ability to play back 
the data in the laboratory for direct examination or further data 
processing. From cost, battery operation and weight considerations, 
the domestic-market portable cassette-tape recorders offer an immed- 
Hope vweatebachiver pOsaibil tev walehousietie cixed 1° 7/8 einch/second 
tape speed limits the frequency response, and tape-transport "wow-and- 
flutter" specifications are far below normal instrumentation-recorder 
standards. 

Again from weight considerations, H-field sensors are not 
attractive, and because of cable weight and installation time, earth 
electrodes for the horizontal E-field components are also unattractive. 
Since the propagating TEM mode leaves the vertical component of the 
E-field with a more useful magnitude, it appears that an E-field 
sensor offers the most practical solution. Again for convenience, the 
commonly used elevated whip-antenna has been avoided, and a more 
readily portable "charge--induction-plate" sensor has been used. (The 


details of this type of sensor are fully covered in Section 3.1.) 
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arn OYSTEM FREQUENCY RESPONSE 


Tiswasenotedein section Jal, 2 that ELF stenals down to at 
least / Hz are of interest because of the Schumann resonances. Below 
about 5 Hz, the system frequency response must usually be rapidly 
attenuated to prevent major baseline disturbances due to local space- 
charge phenomena, geomagnetic signals, and possible artifacts due to 


vibration of antennae and cables. Figure 2.1 shows the general problem: 
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(a) With the lower cut-off frequency at 1 Hz. 


(b) With the lower cut-off frequency at 8 Hz. 


FoLg se aul. Low-Frequency Cut-Off For Baseline Stability 


Large low-frequency disturbances as shown in Fig. 2.1(a) may saturate 
the system amplifiers and lead to signal clipping. Careful choice of 
the low-frequency cut-off can produce the much more useful recording 
Cr mio ee. Cb): 

In the "mid-band" region, interference from 60H2 and its 


harmonics is the obvious concern, although a single high-Q notch filter 
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can usually provide an acceptable solution, (providing that sources 

such as 220V above-ground power lines are several hundred yards away). 
Some trace of 3rd-harmonic noise may then be evident, but in the present 
applications an added 180-Hznotch filter has not been considered 


necessary. 
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(4) With’ the 60-HZ notch filter OUT 


(b) With the 60-Hz notch filter IN 


Pater 22 Notch-Filter Removal of 60-Hz Interference 


Above 60Hz, ELF signals of interest exist up to the nominal 
3-kHz band limit, and of course these higher frequencies may be compon- 
ents of lower-frequency waveforms. In such cases the widest available 
band-pass is always desirable from a signal-resolution standpoint. 
However, to get magnetic-tape recordings of ELF signals below 50Hz, a 
frequency-modulated carrier system is required, and this leads to 
practical limits on the upper modulation frequency. 

As discussed more fully later (under Section 2.4) an F.M. 


carrier frequency of 1000Hz has been used in the system, giving a 
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nominal 200-Hz upper frequency limit for recorded ELF signals. This 
can still provide considerable resolution of the details of a 10-Hz 
transient waveform, so that most of the phenomena related to the 
Schumann resonances are well within the system pass~-band, if 5 Hz and 


200 Hz are taken as the nominal frequency limits. 


te) TRANSIENT RESPONSE 


It was noted in Section 1.1.2 that ELF signals are predomin- 
antly transients generated by far-distant lightning flashes, with 
essentially a random-time distribution; steady-state sinusoidal inputs 
are the exception rather than the rule. 

Unfortunately, for the faithful reproduction of transient 
signals, a limited system frequency response is the very thing that 
should be avoided. The limitations on the time-domain response, for 
various frequency-domain limitations, are well defined by the Laplace 
transforms, but some of the waveform distortions stemming from rather 
complicated transfer functions may not be intuitively obvious, and the 
generation of pass-band artifacts by certain types of input waveforms 
should be specifically noted. 

Figure 2.3 displays some of the implications of a restricted 


low-frequency response and the artifacts generated by typical filters. 


Pulse inputs are shown in the top row (square waves at several frequen- 
cies), the second row shows the output after the signal has passed 
through a typical 2nd-order high-pass Butterworth filter with a 10-Hz 


lower-corner frequency, while the bottom row shows the further result 


of including a 60-Hz passive Twin-T notch filter in the signal path. 
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Fig. 2.3 Typical Signal Artifacts Due to Filter Transients 


(see text for input-output conditions) 


The above illustration of course, merely reminds one that low- 
frequency pulses have high-freauency edge components, that a high-pass 
filter will differentiate the edges of square waves, and that a Twin-T 
notch filter will compound the problem. In both cases, such circuits 
can generate signal "artifacts which may seriously alter the waveforms 
of genuine pass-band signals. 

We might note that the transient responses of standard filter 
types ( Butterworth, Tchebycheff and Bessel ) are well documented [52], 
[53] but in the communications field the steady-state frequency response 


is usually of more concern than such transient responses. For the 
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reception and recording of ELF signals, the reverse is true. 

However, the willingness to tolerate a non-ideal frequency 
response does not really solve the problem. The attenuation at the 
low-frequency end of the pass-band (as noted in Section 2.2) must be 
greater than 20 dB/decade, which implies at least a 2nd-order high- 
pass filter. The most that can be done to minimize time-domain "ringing" 
is to ensure that only real poles exist in the filter transfer function. 
Figure 2.4 shows the limit of improvement by this route; waveform (b) 
shows the output of a 2nd-order Butterworth high-pass filter with a 10-Hz 
corner-frequency, and a l-liz square-wave input signal. The complex poles 
in the transfer function lead to excessive overshoot before baseline 
recovery. Waveform (c) is the output of a 2nd-order high-pass filter, 
with the same corner-frequency of 10Hz but with only real poles in the 
transfer function. The reduction in signal overshoot is evident, but 


the major pass-band "artifacts" of course remain. 


(a) 


(b) 


(c) 


Big 92.4 High-Pass Filter Transients 


(see text for details) 
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The problem associated with the upper-frequency limit relates 
to the system impulse response. For input pulses with a period shorter 
than the filter time constant, the output will simply be some approx- 
imation to the system impulse response, regardless of the fact that 
the pulse "1/T frequency" may be far beyond the pass-band. 

Again, complex poles in the transfer function controlling 
the high-frequency roll-off can lead to overshoot and ringing, while 
maintaining only real poles will minimize the problem. Figure 2.5 
illustrates these points; waveform (b) shows the signal overshoot 
and ringing which constitutes the impulse response of a 2nd-order 
low-pass Butterworth filter, with a 500-Hzupper corner—frequency. 
Waveform (c) shows the improvement when a 2nd-order filter, with 
the same cut-off frequency but only real poles, handles the same 


input pulse, 
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Fig. 2.5 Low-Pass Filter Impulse Response 


(see text for details) 
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The transients introduced by the typical Twin-T notch filter 
can in fact be avoided. This is possible because the 60-Hz interference 
normally is steady state, and the notch filter steady-state response 
can be isolated from signal transients. Appendix A shows, in some detail, 
how a Twin-T notch filter can be modified to give a "clean" transient 
response, while Appendix B offers almost as effective a solution using 
amsoeca lledi. eyratorai notch filter. 

From the foregoing comments, it is apparent that the general 
design philosophy must be to minimize system transients by using the 
lowest acceptable attenuation rates at the frequency-—pass-band limits 
(although total avoidance of complex poles is not really practical) 


and by using a transient-free 60-Hz notch-filter system. 


ZN MAGNETIC-TAPE RECORDING 
2.4.1 DIGITAL vs F.M. RECORDING 


The decision to use a cassette-type portable magnetic-tape 
recorder for field recording of the ELF signals places immediate 
restrictions on the choice of recording methods and formats. 

The magnetic-tape cassettes conform to international standards 
derived from the original Philips Corp. designs [54] which established 
the tape speed at 1 7/8 imches/sec.,( 4.76 cm/sec). This provides an 
acceptable compromise between recording time and frequency response 
for domestic use, but such a low tape speed is an unfortunate restrict- 
ion from a data-recording standpoint. 


Direct analog recording of ELF signals down to 5HzZis not 
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feasible with a magnetic-tape recorder because of the low-frequency 
limitations of the play-back head. The signal output level of a magnetic- 
tape "reproduce" head is proportional to the rate-of-change of the 
magnetic flux cutting the head windings, and this of course is a funct- 
ion of the frequency of the recorded signal, as well as the design of 
the head gap, core permeability and the number of turns of the coil 
winding [55]. The result is that the generated output voltage drops off 
with decreasing frequency at 20 dB/decade. With the available maximum 
tape retentivity of about 1000 gauss and the best available head design, 
the lowest useful recording and playback frequency is about 100 Hz, 
and at 10 Hzthe recorded signal cannot usually be recovered from the 
background noise due to the magnetic-oxide and surface imperfections. 

The use of digital-recording techniques is also ruled out by 
the 1 7/8 inch/sec tape speed, in this case due to high-freauency 
limitations. For digital sampling and a somewhat crude reproduction of 
a sine wave, as few as 5 samples/cycle may be tolerated, but even at 
such a low sampling rate, for a signal at 200 Hzand for an amplitude 
BecOlution ot cay 0.54 (8 binary bits pér sample ) this demands a 
Serialerecoraine bit trate’ of: 

Bit Rate = 9200 x85 x 38 = 78-000) bits/sec. 

In recording digital data on magnetic tape, initially-sharp 
transitions between magnetization directions (corresponding to a 
binary 1 or O ) are blurred by self-demagnetization forces, so that 
the al transition is spread along a finite length of tape. As a 
result, the number of flux reversals per inch of tape which can be 


reliably recovered is not just a matter of reproduce-head gap resol- 
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ution but of tape bit-density capability as well. The theoretical 
limit for the highest-quality magnetic oxides is usually quoted as 
about 5,000 flux reversals per inch [56] but for the data-recovery 
reliability required by a digital system a recording density of 
apouteus000tbits/imchwisva practical@limits"Atrthis bittdensity, 2 
tape speed of 8 inches/sec. would be required for a serial data rate 
of 8,000 bits/second. Since the cassette tape-transport speed is a 
fixed 1 7/8 inch/sec. we would in fact be restricted to an upper 
signal frequency of: 


bbe ps Ba Reote <oLO00 = 7468 Hz 
Sy os vs 


which is about one order of magnitude too low to be acceptable. 

The above limitations therefore restrict the choice of 
recording methods to some form of frequency modulation. An F.M. 
recording technique permits the selection of an audio carrier freq- 
uency well within the record/reproduce range of the recorder to be 
used. The permissible modulation ranges from D.C. to perhaps 204 
of the carrier centre frequency [57]. Such techniques have been 
widely used for many years, and linear voltage-to-frequency converters 
are available commercially or can be designed to specification using 


discrete components. 


2.4.2 PULSE-DURATION MODULATION (PDM) 


The direct use of a frequency-modulated sinusoidal carrier 
is in this case undesirable, since during playback spurious frequency 


shifts due to tape-transport speed variations can be generated, and 


A. 
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such frequency shifts cannot be distinguished from those due to the 
genuine modulation. Tape-transport speed variations of from 1% - 3% 
are typical in portable magnetic-tape recorders ( due to the D.C. motor- 
speed control limitations and fluctuations in winding tension in the 
tape cassette ) and would produce equivalent errors in the demodulated 
analog signal. 

However, one type of pulse-width modulation offers a built-in 
monitor for such speed variations: the fixed-repetition-rate Pulse- 
Duration Modulation (PDM) with "trailing-edge" modulation [58]. This 


modulation method involves the following pulse waveforms: 





| | EAL Se 
Dial? ye Fr aR 
Ty max 
j | 
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(a) Zero Modulation (b) Maximum Modulation 
Fie. )2,.6 Pulse-Duration-Modulation Waveforms 


The "zero-modulation" waveform is shown in Fig. 2.6(a) and has a 50% 
duty cycle ( Ty = 5 ). The repetition rate is held constant, while the 


relative durations of tT, and T, are modulated as shown in Fig. 2.6(b), 


Lees ( Ty + T = =—a Constant 
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and modulation of the transition time of the trailing edge of Ty must 
besiimited to prevent pulse “overlap', that is, qT, or qT, must never be 
reduced to zero duration, or modulation information will disappear, the 
basic repetition rate will be lost, and the demodulator output will be 
meaningless. 

Prewasenoted sins sections 2.4.) that for standard maznetiec tape 
there is a theoretical resolution limit of about 5,000 flux reversals 
per inch. At a tape speed of 1 7/8 inches/sec. this implies a minimum 
pulse width of about 100 microseconds. With a portable tape-recorder 
upper-frequency limit of about 8 kHz,a reliable resolution of 150 micro- 
seconds with an uncertainty of a few microseconds appears to be a 
DEAelLca iaelimit. 

Recovery of the analog signal from the PDM waveform is effected 
by a demodulator whose output tracks the succeeding variations of 


( de aiea ) with the durations of aa and T, measured each "cycle" by 


2 
some form of sample-and-hold integrator. If a 5-sample approximation to 
a sine wave is again assumed to be "acceptable", then a 200 Hz maximum- 
modulation frequency requires 1,000 samples/second, which implies a 


PDM repetition rate of 1,000 pulses/second. This therefore defines the 


pulse widths of the PDM waveform as: 


Ww = i iatlbisieecemel 
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500 usec +350 usec. 


(150 usec min,, 850 wsec max.) 
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Dead RECORD-PLAYBACK OF PDM WAVEFORMS 


Although the pulse-duration modulation offers useful propert- 
ies, the "transmission" medium in this case is a magnetic tape, and 
satisfactory recording and recovery of the PDM waveform remains a 
problem. 

The modulation information of the PDM signal resides in the 
transition time of the trailing edge of iprelative. to its leading edge. 
As already noted, a square-wave signal recorded on tape as an initially 
apbfuptechansesinythe direction of magnetization, does not retain its 


Msquareness"’. 
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Figure 2.7(b) is a representation of how the linear distrib- 
ution of tape magnetization is modified by the self-demagnetization 
forces, so that the original abrupt change imprinted by the record head 
reverts to a relatively-wide "transition zone", whose width is a funct- 
ion of the magnetic-oxide properties and its coating thickness. 

The reproduce-head output voltage during tape playback is shown 
in Fig. 2.7(c), and this signal can then be differentiated to recover a 
signal whose "zero-crossing" points should represent the abrupt time 
oe of the original recorded waveform. In theory, the zero- 
crossing detector could therefore be used to generate a replica of the 
original PDM waveform. In practice, this approach tends to result in : 
poor signal-to-noise ratio, since small amplitude variations of the play- 
back signal and associated tape background noises are greatly accentuated 
By the ditterenttatten, and a zero-crossine “time-jitter" results. In 
addition, for pulse periods much longer than the transition zones, the 
head output voltage has multiple regions of zero slope, which leads to 
ambiguous output signals from the differentiator, hence false zero- 
crossings and output-waveform errors. (The addition of hysteresis and 
some form of voltage-level logic can reduce these false zero crossings, 
but such methods tend to have critically-narrow tolerance margins for 
changing ‘circuit “‘parameters.) 

A much more tolerant system results if the PDM signal is differ- 
entiated pelore it aism@recorded. As) indicated ini Fipure 2.8(c), the tape- 
magnetization pattern will then have a continuously-decreasing intensity 


between each transition zone. The playback-head output voltage will 
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display a zero crossing which is actually related to the end of each 
transition zone, but since no added differentiation is involved, this 
zero crossing will be relatively stable in the presence of associated 
noise. In addition, the slowly-changing intensity of magnetization 
between transition zones provides a residual non-zero output signal of 
sufficient amplitude so that false zero crossings are avoided. 

This method of recording and playback of the PDM waveform 


has therefore been used in the prototype system being described. 


+ 
; 
(a)PDM Signal : | ; | | | | aoe 
(b)Record—Head 
a ee 


Current 

(c)Magnetization 
DivectrLon and Sy es: 2 hae =X 
Intensity icine eee distance 


along tape 


(d)Signal from | ae eee ee 
Playback Head ==Ay/! 
| 
(e) Zero—-Crossing ae 
Detector Output 


Fig. 2.8 Modified Record-Playback System Waveforms 
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2.4.4 COMPENSATION FOR TAPE-SPEED VARIATIONS 


The demodulation of a PDM signal involves the generation of an 


output voltage of the form: 
WA PZ AAG IE) qT, ) 


(with Ty and qT, as in Fig. 2.6) and involves a direct time-measurement 
demodulation. The signal-to-noise ratio of such a system will be related 
almost entirely to the "time-jitter" of the zero crossings of the recov- 
ered PDM waveform. In practice, the separation of the modulation from 
the PDM-carrier waveform is best achieved by a so-called "integrate- 
and-dump" system [59] which is one form of a synchronous detector. 
Although the period iy is normally varied in response to mod- 
ulation, it will also vary due to changing tape speed. However, any 
variation in the total period T must be due only to a tape-speed change. 
Therefore, if the output signal is normalized by a factor proportional 
to the total-period T, it is theoretically possible to make the output 


essentially immune to tape-speed variations. 


lige. 18a is Toe) 15) 


then a fractional change of tape speed will modify the time scale by a 


£ACtOL (0), SO that 


ie = K¢ aT, - af, ) 
aT 
which remains = K( tT, Sel ) 
4b 


and hence no change in output voltage occurs.(The proviso of course is 


that the change in tape speed be slow relative to the sampling rate, 
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So that the same scale factor (a) will apply to Tj, Ty and T. Since 
“wow-and-flutter'" speed variations are usually at less than 10 Fz and 
since the sampling rate is the 1000/sec. pulse-repetition rate of the 


Eitimcarriel tae factor .c) wiliein tact pe essentially constant.) 


Pa) SYSTEM SENSITIVITY AND GAIN REQUIREMENTS 


In a system design, it is of course necessary to have an 
estimate of the required input sensitivity. In the case of ELF signals 
this can be calculated from published data on the measured spectral 
density, of BLM nose. Erom aenumber of sources [60],[61]>[62] the 
average spectral density in the 5-Hzto 200-Hz range can be taken as 
about 0.05(mV/m)2/Hz for the vertical E-field component, but peak or rms 
values enor simply related to the average values derived from 
spectral density data (due to the impulsive character of ELF signals). 
However, for design purposes the expected peak values can be taken as 
perhaps twice the "average" value. 

For a nominal 200-Hz bandwidth, the peak values of background 


noise can therefore be estimated as 


er) aa 0.05(mv/m)2 x 200Hz =» 6 mV/m 


HZ 
with regional lightning storms producing peak signals one or two orders 
of magnitude higher. Using these guidelines, the system requires a "full- 
scale" sensitivity of about 10 my/m on the most sensitive range, and 
gain switching to allow for an attenuation of approximately 40 dB. 


The amplifier-gain requirement is somewhat more difficult to 
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define, since it is dependent on the E-field sensor transfer character- 
istic, and on sensitivity of the pulse-duration-modulator circuit. The 
general requirement is simply that an ELF disturbance of 10 mV/m, whose 
frequency is within the pass band, must provide maximum modulation of 
the PDM-carrier waveform. 

The required amplitude of output signal from the PDM modulator 
PQuLiem tape srecormer tsi typically 1 Volt (peak) inte the “AUX” input. 
However, the fact that the signal is in this case a differentiated 
pulse must be recognized. The automatic-gain-control circuits used in 
many portable recorders are designed to hold a proper recording level 
despite wide variations of the input-signal level, and are of course 
keyed to the average values of typical voice and music waveforms. With 
the low duty-cycle of voltage spikes generated by differentiating the 
PDM square waves, the automatic-level-control setting of the signal- 
amplifier gain can be excessive, leading to waveform clipping and an 
improper recording level. However, if the recorder is also equipped 
with a recording-level meter and a manual level control, satisfactory 
recording is readily accomplished by setting the indicated recording 
level about 10 dB below the level recommended for speech and music. 

On playback, a signal from the "EARPHONE" jack of the magnetic-— 
tape recorder is used to drive the PDM demodulator. The output amplitude 
is adjustable by the recorder "VOLUME" control, with a maximum of 
perhaps 3 Volts (peak) across a typical 10-ohm earphone load. ( Note 
that in some recorder designs, a 10 to 50 ohm resistor may have to be 
provided as a nominal substitute load if the signal is taken from the 


earphone jack for use by high-input-resistance devices. Otherwise, 
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internal negative-feedback circuits may display instability, so that 
recorder-amplifier oscillations may occur.) 

The demodulator output voltage is a matter of choice and 
design. In the system being described, the output has been scaled to 
45.0 volts full scale, for compatability with an available analog/ 


digital computer input. 
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CHAPTER 3 


DETAILED SYSTEM DESIGN 


The general system configuration is shown below, in "block- 


diagram" form: 
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(b) Playback Configuration 


Figure 3.1 System Block Diagram 
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The detailed design equations, circuit schematics and notes 
on circuit operation are included in the following sections. The discuss-— 
ion can be conveniently subdivided into four functional sections: the 
E-field sensor assembly, the signal amplifiers and filters, the PDM 
modulator, and the PDM demodulator. 

It will be noted that recordings of system waveforms have been 
extensively used for illustrations. The ease with which such recordings 
could be made, with almost any desired time scaling, resulted from the 
fortunate availability of a Hewlett-Packard Model 5480A Signal Analyzer. 
This unit provides analog/digital sampling, with data memory-storage 
and playback facilities. A "recording™” can be made consisting of 1,000 
sequential samples at a choice of sampling rates from 2 samples/second 
to 100,000 samples/second, and "played back" at a much lower rate via 


digital/analog output facilities to a standard X-Y recorder. 


GL ThE E=P PEED SENSOR 


For system Seeisvueten the usual E-field sensors mentioned in 
Section 1.2.2 are considered rather inconvenient. The following sub~sect- 
ions cover the design of a "charge-induction" type of sensor, which app- 
ears to be a useful alternative to the whip antenna for measuring the 


vertical E-field component of ELF signals. 
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3.1.1 THE CHARGE-INDUCTION PLATE 


One of the simplest methods of measuring the vertical E-field 
at the Earth's surface is a "charge-induction-plate" method, which is 
basically the system used by Wilson in 1916 [63] to measure the Earth's 


surface-charge density. 
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(a) The Wilson "Charge-Plate" (b) The Portable E-Field Sensor 


Figure 3.3 The Charge-Induction. Plate System 


Wilson took great care to keep the charge plate flush with the 
Earth's surface, and used a null-balance method to prevent E-field 
distortion (since he was concerned with absolute charge-density measure- 
Ments )weforea portablersystem this is impractical, but as in Fig. 3.3(b) 
above, a shallow above-ground assembly has been used. This leads of 
course to some edge distortion of the incident E-field, but since the 
system can be calibrated in known fields, this distortion factor is of 
little concern. Note that the "guard" electrode is in this case used 
to reduce the effective mounting capacitance, not to eliminate electrode 


edge effects. (See Section 3.5 for the physical design details.) 
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3.1.2 THE SENSOR EQUATIONS 


The signal voltage van available from the charge-induction 
plate for any incident E-field can be calculated from basic charge- 


capacity-voltage relationships. 
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Figure 3.4 Equivalent Circuits For The E-Field Sensor 


In Fig. 3.4(a) the normal incident E-field induces a surface-charge 


density of 


iC. is the free-space capacity of the exposed plate, and Ci, is the 


lumped capacity of the plate mounting, wiring, and electrometer input, 


then 
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vo Gee Ca 
Gn ft Cin 
This instantaneous voltage will decay exponentially as the induced 
charge leaks off through associated leakage resistances represented in 
Fig. 3.4(b) by Rj,,. 

For system analysis, the equivalent circuit of Fig. 3.4(b) is 
convenient. The incident vertical E-field can be viewed as the result of 
aavoleaze cource Vavappiied to a theoretical) sheet-electrode located one 
metre above the sensor plate. The capacitance C, remains valid and can 
be calculated directly as a parallel-plate condenser, and any E-field 
can be defined directly in volts/metre. 

Using "standard" Laplace notation, the voltage transfer 


function can be written as 


ee) Ca \ Ss 
H(s) We) Ges Cin/ s + ae 
Rin(Cg + Cin 


The charge-plate system is therefore a lst-order high-pass filter, with 
a lower corner-frequency of 


if 
u) SS 
in cer ic) 


and for pass-band frequencies (for w>>wW .), the electrode system will 


appear as a simple capacitive voltage divider with 


EOL ea eee ne oe 
Va Ca + Cin 


For the 10 inch by 6 inch charge-induction plate used, the 


value of C, is approximately 0.33 picofarads, and the total anput 
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Capacitance of a typical electrometer circuit can be taken as approx- 
imately 10 picofarads. Therefore, for an E-field change of say 10 mV/m, 


the voltage signal available at the electrometer input will be 


G We 
Vee x Vv a = Vv fee oe : Vv. 


( Greater sensitivity could of course be achieved by increasing the 
dimensions of the charge-induction plate, which would increase the 
value of C,.) 

If desired, the amplitude-response versus frequency can be 
calculated from 


7 a C eee ctien la ae OS 
lH(jw)] = Gw) ee) | 1 + (wo/w)> | 


L 


where UD a a ge aaa coo ata: a Vim an 
Ca C in) 


The transient response to a step input (such as the leading or trailing 


edge of a square wave) can be given in the time domain as 
C 
ee ee ae oe fee P 
v(t) Co; Z a exp[-t/Ri, (Cz + Cin)! 


For a lower corner-frequency of say 1Hz it is evident that R,;_ 


must be a very-high resistance, i.e., 


Re. = s = AAS 10!° ohms 


pees ee eee 
in Pomp CAe ace) be Lenw( 241) GL0= "> 


which is greater than 10,000 megohms, but is a Fealdzablemmput resict— 


ance provided by many alternative electrometer circuits. 
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Sees THE ELECTROMETER-AMPLIFIER CIRCUITS 


As noted in Section 3.1:2, for the desired low-frequency 
response and the necessary voltage sensitivity of the E-field sensor, 


the following input conditions appear necessary: 


(Ann G. el Opileorarad 
eae minimize 
10 
Ra el) ohms 


The simplified schematic of the electrometer-input circuit is shown 
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Figure 3.5 The Electrometer-Input Circuit 


This configuration is essentially a unity-gain source follower, with a 
"bootstrapped" input resistor R,. Since large electrostatic voltages 


can be encountered during handling and use of this "charge-induction" 
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type of sensor, gate protection of the input FET (field-effect transist- 
or) is imperative. Input capacitance must be as low as possible, and the 
1/£ noise figure below 10 Hz is of concern. These considerations and 
component availability have led to a number of design compromises, alth- 
ough the electrometer performance remains adequate. 

Transistor Q, is a dual-gate MOSFET (an RCA 40822, which has 
integral gate-protection diodes). This device retains a useful transcon- 
ductance at the very—-low drain current (100A) used to minimize its 
noise figure, and the second gate provides isolation between the input 
gate and the drain, to reduce the effective input capacitance. Transistor 
Q5 is a low-noise type 2N5087 PNP device, also operated at a low collector 
current (100uA) for reduced 1/f noise. 


For such a circuit, the voltage gain is approximately 


ee Emi beerRy 


f 
Te misao" 4 


and is therefore maximum for the maximum Ry product. For the type 


Sm] 
40822 MOSFET with zero voltage between gate 2 and the source, and with 
Vas = 3 V, the best compromise between noise and gain appeared to be 
with ty = 100uA and Ry, = 5K ohms. Under these Conditions 2.74 = 107"mhos 
and Rw iti sie LOOWA the 2N5037 has an he, = LOOMPsoOuLnat actiad cain 
of the "unity-gain" amplifier is approximately A, = O226s 

The RjC, time constant of the “bootstrap” coupling network is 
kept greater than 1 second to ensure its low-frequency effectiveness to 


at least one decade below the nominal 5 Hz low-frequency cut-off of the 


over-all system response. 
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Themiaptverestetor ky ts a Victoreen Corp., Hi-Meg” electrom- 
eter resistor of 101° ohms. With the "bootstrapping" arrangement, the 


effective input resistance will be 


which is more than adequate in this application. 


these leckbRomeberssec &1OnLOm.F1e. 03.5 2S A.G..,.coup led, toa 
cable-iyive amplitien, witch is sketched in Fic. 3.64 [hie provides a 
nominal voltage gain of five, and provides a low-impedance cable drive 


to minimize noise pick up. 
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Piguress. 6 The Cable-Drive Amplifier. 


Transistors Q3, Ons Qn are all low-noise units, type 2N5210. For minimum 


1/f noise at 10 Hz, Q3 is operated at a collector current of 20 uA, and 
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and the optimum source resistance of 5K ohms (for minimum noise at this 
current level) is provided by the resistance of R5 and Rg in parallel. 
The emitter resistor Ro provides bias stability, and is left unbypassed 
for flat low-frequency response. 

Withesuchya slow collector current, the output resistance of Q3 
ds) tairiy hich. and the collector—load resistor R7 must be high ( 120K 
ohms) to achieve an acceptable gain. Transistors 01 and Qs constitute a 
Darlington-connected emitter follower, which is necessary to minimize 
the loading on Q3 and which provides the low source impedance for 
driving the shielded cable. Effective output impedance is approximately 
10 ohms. 


Capacitor C3 provides an initial by-passing of radio-frequency 


and high audio-frequency noise, with the roll-off starting at about 2 kHz. 


A direct-current feedback loop from output to input is provided 
by Ry, to provide low-frequency stability, and sets the net voltage 
gain, from sensor plate to cable-drive amplifier output, at a nominal 
ValuerOr ..8 (ltaimont beinoted that with hicher gain, the single= 
shielded cable permits sufficient signal leakage back to the very-high- 
impedance sensor plate, so that low-frequency oscillations can occur 
in the 1-Hzto 10-Hzrange. Double-shielded cable would presumably elimin- 


ate this problem.) 


The detailed circuit schematic with component values used, 


is given as Figure 3./. 
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5.144 SENSOR CHARACTERISTICS 


The characterization of the sensor assembly requires the meas- 


urement of Rap and Ci, This is effected by applying a square wave of 


known amplitude first to the sensor plate directly, and then to a sheet 


conductor supported a known distance above the sensor plate. The output- 


signal amplitude and waveform can then be used to calculate ee aan 
and the net gain from sensor plate to the output of the cable-drive 


amplifier. The guard-electrode extensions prevent distortion of the 


incident ™H-field. and ane at D.C. “eround potential. 


sheet conductor 


SS ss 











guard ~ ae Sah 

extension ee 
cable to 
amplifier 


Figure 3.8 Sensor Calibration Arrangement 


The results from such a series of tests are shown in Figure 3.9, 


The input 5-Hz square wave with an amplitude of 25 mV (p-v) is shown in 
Pipes .o(ale LE this ssdenal is) connected: directly to the sensor plate, 
the cable-drive amplifier output is as shown in Fig. 3.9(b), from which 


it is evident that the sensor assembly net voltage gain is 
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(a) V, (20 mV/cm) 
(b) ee ae ee V, (200 mV/cm) 
(c) Vo (50 mV/cm) 


(x = 50 msec/cm) 


Figure 3.9 Sensor Calibration Waveforms (see text) 


With a sheet conductor mounted 10 cm. above the sensor plate, 
the capacity C, (om the vequivalent circuit 3.4(b)) can be calculated 
assuming parallel sensor-plate areas, 10 cm. apart, giving C, = 3.3 pF. 
Tf the input sienal is connected to this overhead sheet as Vee the 
measured output of the cable-drive amplifier is as shown in Fig. 3.9(c), 
from which Vo = 55 mV. From this, knowing the amplifier gain to be X 8 


i eee On Ce 





Then, since Vimy = eel. 


Mee 
and rearranging, Cin = Cag — Vin) 


then Cin 


(Seo tot) == aeRO eompias 
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After the step response, the waveform of Fig. 3.9(c) shows 
an RC discharge curve due to the combined effects of the sensor time 
constant and the amplifier coupling-capacitor Co Cterie.. 2.6, From 
the difference between the indicated time constants of Figures 3.9 
(ip eamca CC.) at ie exponential decay is estimated to be 0.44 of the 


step amplitude in 100 msec. 





Then, since e 0 L/RC = 0.44 
RC = es =20) 1 2esec. 
which is Ree Game C ines L2usec. 
giving Rin = Fee = 101° ohms. 


(It might be noted that although the effective resistance of the 1019 
ohm input resistor should be about 5 x 10!! ohms due to the bootstrap 
circuit, this will be shunted by the leakage resistance of the gate- 
Proeceeion diodes in the 40822 MOSPET. The effective Ri, = 101° ohms 
as calculated above, is therefore the resultant of these paralleled 


resistances.) 


A noise figure for the sensor input system can be derived from 


the following test records: 


IGN Ge O ohms 





(Hor. sceale 
= 100 ms/cm) 


(b) PEN PON otf arn gle A a Pago get nn Bu et = 1Otachine 


Figure 3.10 Equivalent-Input Noise Records (see text) 
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Figure 3.10(a) Bhawe the amplifier-output signal when the sensor plate 
is shorted to ground. The peak-to-peak amplitude of this noise is approx- 
imately 30 wV referred to the input (or approximately 7 uV rms for 
gaussian noise, assuming the 20 level as the nominal peak value). With 
the electrode non-shorted but electrostatically shielded, the noise 
level is as shown in Figure 3.10(b), and is approximately 80 uV (p-p) 
or 20 pV rms, referred to the input. 

The theoretical noise level can be calculated on the assumption 
that the dominant source will be the Johnson (thermal) noise in the 


101° ohm input resistance. The classical expression for Johnson noise 


anes 
Vi, = 2\kIRB Volts (rms) k = Boltzmann's constant 
T = temperature, K 
R = resistance in ohms 
B = bandwidth in Hz 


In this case, the noise bandwidth is limited by the Thévenin-equivalent 


Caigeuats 





Figure 3.11 Thévenin-Equivalent Circuit for Input Noise 


which sets an effective upper corner-frequency at approximately: 


LS : 


oo 18 
CT Riatoe cia) 2nC10 i 0-82) ; 


The slower cornér-prequency for the noise is at 053 Hz, set by the RG 


coupling network into transistor Q3 of the cable-drive amplifier. The 
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effective noise bandwidth is therefore a nominal 1 Hz, and the expected 


thermal-noise level will be: 





<4 
i} 


Pee Ceo 07 2808 x 102) (1077) (1) 


i 


palates 
2A Wee tema 13 x 1055 213 uv (rms) 


which is in reasonable agreement with the observed noise level of 
about 20 uV (rms). 
A sensor noise figure can therefore be defined and calculated 


as: 


nee 10 1 total noise-power output 
: 08 noise-power output due to Ry, 


i} 


\ 2 
WO Tae fa 5 33) GP 


which indicates that the sensor noise level will be 3.8 dB greater 
than the theoretical noise level due to thermal noise in the very-high 


input resistance. 
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Se AMPLIFIERS AND FILTERS 


As noted in the system block diagram of Figure 3.1(a), part 
of the field-recording system consists of the signal amplifiers, a 60-Hz 
notch filter and a high-pass filter with a lower corner-frequency of 
about 5 Hz. These are arranged as sketched in Figure 3.12 below, and 


will be discussed in detail in the following sections. 


















Signal High-Pass 
from Amplifier Filter to PDff 
E-Field Modulator 
Sensor 









Attenuato 60-Hz 


Amplifier 


60-Hz Notch Filter 


Figure 3.12 Amplifiers and Filters: Block Diagram 
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3.2.1 THE INPUT AND SUMMING AMPLIFIERS 


The amplifiers are designed around integrated-circuit operat- 


ional amplifiers. The simplified schematic is shown below: 
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Figure @2i3— The Input ‘and ‘Summing Amplifiers 


The effective input-coupling time constant RC, is selected so 
that its low-frequency roll-off starts at about 1 Hz, considerably below 
the 5Hzdominant lower corner-~frequency set by the high-pass filter 
(Section 3.2.3). Resistor Rj is a relatively-low 22K ohms, so that the 
input-resistance changes of amplifier Al due to attenuator switching 
will not materially alter the input time constant. 

The gain switching of amplifier Al is defined in terms of 
attenuation for system-application convenience. The resistance values 
of Ro,R3,R,,R5,Rg are selected so that voltage-gain changes are in 


10 dB steps, as tabulated on the following page: 
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Attenuator Voltage 
Setting | _Gain 
O dB X 100 
=1') <dB He Seb 
= 0 aD x 10 
—s08 db Ker 
-40 dB Joe 


The voltage gain of amplifier A2 is set at X 2.5 by the resist- 
ance-ratio Rg/R7. This amplifier also serves as a "summing" amplifier 
for the normal signal voltage Vj and the out-of-phase 60H2 signal Vo, so 
that the 60-kz interference is effectively nulled out. 

Capacitor Cy across the feedback resistor of amplifier A2 
provides the dominant high-frequency attenuation, with the roll~off 
starting at about 200Hz.A small D.C. feedback is provided by resistor 
Rg to minimize output D.C. shifts with different attenuator-switch 


settings. 


O22. 2ee ee OU-HZ NOLCHAELLIER 


It has already been noted that 60-Hz power-line interference is 
within the ELF pass band of interest ( 5Hz- 200Hz). Although system 
portability does offer the opportunity of selecting a recording site 
far removed from overhead power lines, the fact remains that 60-Hz 
interference can be so all pervasive that even in such locations a 60-Hz 
notch filter may still be required. 

It was specifically mentioned in Section 2.3 that the transient 
response of the typical Twin-T notch filter can generate signal artifacts 


if the input waveform approaches the shape of a square-edged pulse. The 
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following notes offer amore detailed discussion of this problem. 

The "standard" passive Twin-T RC notch filter is shown in 
Figure 3.14(a), in a symmetrical configuration. Figure 3.14(b) shows 
the filter transient response for various input pulse widths, for a 


notch frequency of 60 Hz. 


Input Pulse 
Width 


ARAN 100 msec. 


a eee 10 msec. 


| iamsecr 
Givaieaescry 








Schematic Transient» Response (VG.,_) 


(a) (b) 


Figure 3.14 The Passive Twin-T RC Notch Filter 


With the component symmetry noted in Fig. 3.14(a) and with no 


output loading, the voltage transfer function is 





MOU ome Set Wyo 
See FS H Ss = 2 Bw Die ya) 


wheres we = L/(RG). 
From the above transfer function, it is evident that the Q of 
this configuration is only 0.25, so that for a 60-Hz notch frequency, the 


rejection-band width is 240Hz ( putting the -3 dB points at 16Hz and 
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22582). This notch width “is of colrse unacceptable when the desired pass— 
band*is 5 Hz to 200%z, and the signal artifacts of 3.14(b) are certainly 
undesirable. 

Active network high-Q versions of the Twin-T notch filter are 
well known [64],[65] and can reduce the notch width to almost any desired 
value. However, the poor transient response remains and high-O ringing 
can compound the problem. 

The poor transient response as illustrated in Fig. 3.14(b) is 
related to the fact that the desirable notch in the frequency-response 
curve is a steady-state characteristic, while the transient response in 
the time domain is dominated by the RC-coupling configurations and their 
time constants. The input network ( C and R/2 ) simply differentiates 
the edges of applied square waves or pulses. 

Fortunately, in the present application, the steady-state 
response can be almost totally isolated from the transient response. In 
the modified notch-filter arrangement of Figure 3.15, the steady-state 


60-Hz interference is isolated via a sharply peaked 60-Hz amplifier (the 
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high-Q notch filter is now located in the negative-—feedback loop of an 
operational amplifier). This out-of-phase 60-Hz signal is then added to 
the original signal with its amplitude adjusted to "null out" the 60-Hz 
interference. The transient signals pass directly through, and never 
enter the twin-T network. 

With typical component values as detailed in schematic Figure 3.21, 


this modified notch filter has the voltage transfer function: 


ei (=0045 : 045k) , oe (2:01 2 175k) 42 
CS poe (Wee oles) 113i. 2 = See ° 


6-8 0045a0se4+ 4 Ole 


where k, is the fractional setting of the null-adjust potentiometer. It 
is evident that the second term of the bracketed numerator can be 


weducedm@toy zero by adjusting the value of ky ( to about 0.1 ), leaving 


(EIN jieINSIESIS itGMAVeETOim, IG}S 


2 2 
7 Sa 005s 
BOS (x + .0045w,s + at) ES) 


The steady-state zero has been reestablished, although a slight 
shift of frequency is necessary to restore the notch to precisely 60 Hz. 
(In practice, both the frequency and the null adjustements are made while 
observing the filter output with an oscilloscope, with a 60-Hz input 
signal derived from the power line.) The value of .0045 in the denomina- 
Peete weenie reranSter snunrceLOn inaicates an etmective OLlor 1/0045 2220, 
Wiikchesives a, noOtehewidthot 60/0 =.4.2/ Hz at the. -.3 dk points. 

The modified notch filter has a pulse response controlled by 


the characteristics of the summing amplifier, not by the twin-T network. 
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Only a residue of the high-O ringing of the Twin-T network remains 
visible in the over all modified notch-filter response. The pulse 
response for the modified system is shown in Figure 3.16 below, for 


different input~pulse widths. 


sinput-Pulse 


| | “Width. 


100 msec. 
10 msec. 

eek... 
—~—- lsmsece 


Figure 3.16 Pulse Response of the Modified Notch Filter 


A much more complete discussion of the modified notch filter, its 
derivation, and the manipulation of the transfer-function poles, is 
given in Appendix A. 

ttwshouldsbe movedsthat also—-called Weyrator Wnotch—-11i1ter 
circuit is often included in manufacturer's integrated-circuit applicat- 
tonsnotes [66], [6/ |) with no comment on its Special characteristics. Lt 
dpeseim tact have a transter function simidar to Eq'n (3-2.2.b) with 2 
pulse response equal to that of Fig. 3.16. For medium and low-Q applicat- 
ions (Q<20) the gyrator would be preferred to the modified Twin-T 


system because of its circuit simplicity. For high-O applications (0>20) 





ii wath “a, 9 yy 5 





me 


vas 


“ma hee 
CG 1 Peps {5 2 vv 
f| ve 
ey 
f 
7 om En b74§ 
? 
oJ 
| 
My ' 7 
i j 
uit > - i = 
r! hi Wd =@ 


el) 





7 i 


2. wala — 


i 


tha 2d 


- 
’ 


hee | é 


ie toed i 


fd kw tae 








i nan 8 


a a 
ase” ; 


Peal 


undesirable capacitor ratios become necessary, which ruled out its use 
in the present ELF-receiver design. However, for completeness, a brief 


discussion of the gyrator notch filter is given in Appendix B. 


Sey is: THE HIGH-PASS FILTER 


A basic second-order high-pass active filter is shown in 


Fieuress o7(a): 






|Basic High- 
az Pase iPaLlleeie 


Vin 





Basic High-Pass Filter Modified High-Pass Filter 


(a) (b) 


Figure 3.17 The High-Pass Filter Modification 


The voltage transfer function of the basic high-pass filter is: 


Vout (s) 


—_——_—__.. Bg 
Gates ee = 


sg? + [wa + 1/ (C1R9) + (1 - K)w1]s + (W109) 


where wy = 1/(R1C]) and wo = 1/(RgCo). 
To retain only real poles in Hj(s) (as suggested in Section 2.3 
to reduce transient overshoot) component values can be chosen so that 


Rip = Royand Cy = Coeeinen tor an amplifier gain factor-of K = +1, the 
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transfer function would reduce to: 


sg? 32 


Hailes 2 Ue AIRES Rheem (aeten)y2 


where Wo. = (Ra Cap = 1/(R2C2), and the poles would both be real at 

Ss = -Wo. This would provide a low-frequency roll-off of 40 dB/decade, 
but for f, = 5Hz the response would only be down by about 28 dB at 1Hz, 
which in practice is insufficient to maintain an acceptable base-line 
stability in the face of large E-field disturbances encountered at about 
this frequency. 

The best approach appears to be the addition of a transmission 
zero around 1Hz, which can be done effectively by the filter modification 
shown in Figure 3.17(b). For the desired frequency-response curve, comp- 
lex poles must be tolerated ( Ry = Ro/20, Cy = 20Cy ) but the zero 
locations do provide some overshoot damping. (See Fig. 3.18 and Fig.3.20.) 

The transfer function of the paralled low-pass network ( R3 and 


Capote ieure 317 (bh) i ici: 


ee a 
Ha (s) Sis W3 
where w3 = 1/(R3C3) = wo. For the component values used (noted in Figure 


SoZ. yoand with they (x 0nl)) scale factor provided by appropriate resis-= 
tance values of the subsequent summing network, the resultant transfer 
function becomes: 
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or in factored form, 
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These pole and zero locations can be graphically displayed as in Figure 
3.18(b) to permit some visualization of their effect on the system 
frequency response. 
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a ra + 
Basic High-Pass Section Pole-Zero Locations of the 
With Complex Poles Modified High-Pass Filter 
(a) (b) 


Figure 3.18 High-Pass Filter Pole-Zero Modifications 


One penalty for securing the added attenuation rate just below 
5Hein this manner, is the fact that a D.C. term has been added to the 
transfer function ,[0.lw,/(s + w.)] which for a notch depth of about 
-60 dB constitutes a zero-frequency leakage at the -20 dB level. (See 
Figure 3.19.) 

Actually, this can be tolerated since the input coupling 


network provides a 20 dB/decade low-frequency roll-off, and the sensor 
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amplifier adds another 20 dB/decade attenuation from almost the same 


point. The effective low-frequency roll-off is indicated by the dotted 


curve in Figure 3.19. 








modified 
high-pass 
response 


effective 
owere ailliL 
response 








Figure 3.19 Frequency Response of the High-Pass Filter 


iiewattenuation wate Lrom=o-ezdown to 1.5 Hz2as"about 100 dB/decade, and 


provides an acceptable base-line stability, (as was indicated in Figure 


Hs Alb ye 
The transient response is of course somewhat degraded from that 


possible with only real poles, and is noted in Fig. 3.20 below. 
Input Square Wave 
(isHze) 


Output of Basic 
High-Pass Filter 
(real poles only) 


Output of Modified 
High—Pass Filter 


Figure 3.20 Transient Response of the Modified High-Pass Filter 
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Sip) THE PULSE-DURATION MODULATOR 


The ouput of the modified high-pass filter network of Fig. 3.21 
is the input signal to the pulse-duration modulator circuit. The opera- 
tion of the modulator section can be conveniently discussed under a 


number of sub-sections, as outlined below in block-diagram form. 


‘Late Ramp | 















Zero Modulation | Generator | 
Palse wicthesaiust PDM Signal 
/ 
T 
Signal Amplifier Voltage Different- 
Input /Limiter | “(Comparator fOLOT 
|& Scaling | 


oe 





{ 
Wniareieaies iModulation | 
| Indicator Output to 
7OPUE | Magnetic-Tape 
Recorder 


Fiourero. 22 Pulse-Duration Modulator Block Diagram 


The following sub-sections cover design details, schematics, and 


circuit operation. 
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Sha7 ll BUFFER-LIMITER AMPLIFIER 


A buffer-limiter amplifier is provided between the modified 
high-pass filter system (of Section 3.2.3) and the pulse-duration 
modulator. This amplifier performs several functions, and its general 


configuration is shown in Figure 3.23 below. 






Signals From R6 

The High- Ry 

Pass So ANS ONT OSTA 

Filter Ro | To Voltage 
Section Re NANT ee te 8 Comparator 
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Switch 


Figure 3.23 The Buffer-Limiter Amplifier 


The signals from the low-pass shunt network and from the high-— 
pass filter section are combined via scaling resistors Rj and Rp», and 
become the signal input to the non-inverting amplifier A/. 

A voltage gain of about X 13 is also provided by this amplifier 
(set by the resistance ratio (R5 + Rg)/R5). This final voltage gain is 
chosen so that an E-field disturbance of 10 mV/m at the sensor plate, 


with an attenuator-switch setting of 0 dB, will provide a voltage swing 
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of approximately 6 Volts at the output of amplifier A7. This constitutes 
a "full-scale" signal, which will be used to drive the pulse-duration 
modulator. 

As noted in Section 2.4.2, a pulse-duration modulator should 
not be over modulated, or the basic pulse-repetition rate may be lost 
by "pulse overlap". Any danger of such over modulation is readily 
prevented by limiting the signal voltage at the ouput of amplifier A7, 
and since this limiting is non-critical, (it implies off-scale data), 
it is convenient to simply use the operational-amplifier saturation 
limits as an intrinsic signal-limiting system. 

A panel MARKER switch is provided so that an "end-of-record" 
signal can be added to the data being recorded. This marker signal has 
been chosen as an off-scale signal (via Ry) of sufficient magnitude 
to -Over-ride any normal signal. 

ASEnOLedmi ie sect1 on a.2).0,8traces of the D.C. ofEset voltages 
of operational amplifiers Al and A2 will be passed through the low-pass 
shunting network of the modified high-pass filter system. In addition, 
the positive and negative clipping levels of amplifier A7 will differ 
by perhaps one-half volt. These combined offsets are trimmed out by a 
fixed-bias resistor R3 during assembly and calibration. (This is effected 
by observing the output of amplifier A/7 with an oscilloscope, while 
introducing an adjustable sine-wave signal into the system. The signal 
amplitude is adjusted so that a slight signal over-drive is observed at 
the output oLvamplitier Ay, and resistor Be is. then selected so that 


symmetrical peak clipping occurs.) 
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B23 32 THE 1-kHz RAMP GENERATOR 


The pulse-duration modulation scheme involves the use of a 
voltage "ramp", i.e., a voltage which rises linearly with time between 
predetermined) lamits. The voltage is "reset" to its een, veils at the 
end of seach millisecond, so that wits repetition rate is 1 kHz. The 


simplified schematic is shown below. 








= ' 
1-kHz Ramp, 
to Voltage 
Comparator 


Figure 3.24 The L-kHz Ramp Generator 


Amplifier All is arranged as an integrator, with its output 


voltage being related to time by the expression: 


Teansistor Qj sis formally conducting, So Thats the @riective: value of 
Vj, 18 normally constant at the -6.5V saturation voltage of the collector 


OnOn. thie initial (Keset) value of Vo at time t°= 05 16 approximately 


-4.9V, and using the component values noted in the detailed schematic 
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of Figure 3.29, the above equation can be evaluated as: 


(-6.5V)t 
Vo(t) = - aa Os os, | - (4.9v) satceO 


I 


(10*t - 4.9) Volts 


from which it is evident that V,(t) is a positive-going voltage ramp, 
Pic inesaewancoustant rate: of .0V) msec, and starting from an initial 
value of -4.9 Volts. 

During this rise time, amplifier A10 is saturated with its 
output at approximately +6.5V. Resistors Ry and Ry scale this down to 
+4.7V which is then used as the reference voltage on the (+) input of 
AOpechce—)minputceoLeAlOets the ramp voltage V, (rt), monitored via 
resistor R7. 

The ramp is terminated when V.(t) exceeds +4.7V, i.e, when 
the (-) input of amplifier Al0 becomes positive with respect to the (+) 
input. When this occurs, the output of AlN switches abruptly from its 


positive-saturation level tothe negative-saturation level of about -6./7V, 


(which switches the (+) input-reference value to -4.9V). The ramp durat- 
ion is therefore the time required for Vo to rise from -4.9V to +4.7V 


at the rate of 10V/msec, or 


Q 
ab AA 44D) nsec = 960 microseconds 


When the output of amplifier A10 switches to its negative level, 
the base current for transistor O7"drops*to zero,vandeits collector is 


no longer clamped to its -6.5V level. The RC time constant of the All 
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integrator is now much reduced, since the diode Dj conducts, and the 
effective Vj, to the integrator is 46.8V through a resistance of only 
17K ohms, (Ry, at Re). This generates a negative-going ramp, starting from an 
upper-switching level of +4.7V, with the relevant expression for V.(t) 


becoming: 


2 aR (46.8V)t \ 
ot) > = = aa (4.7) Dest 


(-1.8x10°t + 4.7) Volts 


Since this ramp voltage goes negative at the rate of 180V/msec, it 
reaches the -4.9V level (of the (+) reference input of the voltage- 
comparator A110) in the time 


180V/msec = he usieie: 


Ag) 


This constitutes the ramp "reset", after which amplifier A10 switches 


to its original state and the next voltage ramp begins. 










Effective Vian 
to the 
Integrator 


V et) 


Figure 3525 The 1l-kHz Ramp-Voltage Waveform 
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Ent en 8) THE VOLTAGE COMPARATOR 


The required pulse-duration-modulated waveform is generated by 
the ouput switching of amplifier A9, which is arranged as a voltage 
comparator. The ouput-pulse duration is determined by the time taken 
for the dinear ramp to rise from its initial value to the level of the 


modulating signal. 


Voltage 
Comparator 
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Bete ns e220 The Voltage Comparator 


The signal-input voltage (from the buffer-limiter amplifier) 
is scaled down to plus-or-minus 3 Volts maximum, by resistors Rj and Rg 
so that all voltage comparisons in amplifier A9 occur well within the 
linear common-mode range. The l-kHz ramp signal from integrator All is 
also scaled down, by resistors Rg and Ry tO 4. 4V land) 3.6) jamits. 
This ensures that the ramp magnitude is always greater than the modul- 
ation signal, so that the basic pulse can never be reduced to zero 


duration. The ramp voltage therefore sweeps the (-) input of the voltage- 
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comparator from -3.6V up through the signal-voltage level present on 
the (+) input. When the ramp voltage exceeds the signal voltage, the (-) 
input becomes positive relative to the (+) input, and the output of 
amplifier A9 switches abruptly from its positive-saturation level to the 
negative-saturation level. The output reverts to the positive level 
when the ramp is "reset", and the (-) input of A9 reverts to the initial 
-3.6V level. The duration of the positive-output pulse is therefore 
linearly controllable by the signal level at the (+) input of voltage- 
comparator A9. 

The time relationships of the various waveforms involved are 


sketched in Figure 3.27. 














(+) Input 
of AQ a 





—» <—150Us 


Input Signal: Max. Neg. Inputeoienge i: Max s fos. 
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Figure 3.27 Pulse-Duration-Modulator Waveforms 


It has already been noted that the positive and the negative 


saturation levels of the operational amplifiers are slightly different, 
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so that both the input signal (when adjusted for symmetrical clipping 
in the limiter) and the 1-kHz ramp voltage have slight D.C. components. 
The PULSE-WIDTH-ADJUST potentiometer P, (of Fig. 3.26) adds a selected 
bias voltage to the (+) input of the voltage-comparator A9, to control 
the "zero-modulation" switching time. This potentiometer is simply 


adjusted so that Ty = To with zero-input signal, as observed with an 


oscilloscope at the output of A9 or at the output to the tape recorder. 


3.3.4 DIFFERENTIATION AND SCALING 


The pulse waveform at the output of thevoltage-comparator A9 
is a true PDM waveform with trailing-edge modulation, but as noted in 
Section 2.4.3, an improved magnetic-tape record-reproduce characteristic 
is obtained if the pulse waveform is differentiated before use. The 
PDM waveform from the voltage comparator is therefore passed through 
an RC coupling network consisting of Cj and (Rs FeRe), wich basa very 
short time constant (relative to the minimum pulse width of Ty or To). 
The resultant spike exponentials, representing the transition times of 
the PDM waveform, are scaled by Rs and Rg to about +1V peak, and are 
then fed to the output jack, ready for introduction into the typical 


"AUX'' input of a cassette tape recorder. 


3.3.5 THE MODULATION INDICATOR 


The magnitude of the natural ELF signal varies with time, and 


in some cases with the recording site. Some indication of peak field- 
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strength is therefore required during the recording process, so that the 
best attenuator-switch setting can be selected for a reasonable PDM- 


modulation level. 


To Modulator RQ 





Figure 3.28 The Modulation-Level Indicator 


Referring to the above schematic section, note that the signal being fed 
to the PDM modulator is "sampled" by diode Dj], with the positive peak 
values stored on €;. The voltage on Cy is continuously monitored by the 
high-input-impedance voltage follower A8, with a slow discharge rate 
via the input resistance of A8. The +3V maximum signal level to the 
modulator is restored to +6V maximum by a gain of two in the voltage- 
follower (set by the resistance ratio (Rj + R2)/R 1). The low-impedance 
output, via a calibration resistor R3, drives a small panel-meter M, so 
that a full-scale E-field transient produces a full-scale meter deflect- 
LOT « 

This indication of the modulation level is rather "crude", due 
to the random transient nature of the ELF signals. However, it has been 
observed that typical signals have reasonably equal positive and 
negative peaks (on the average), so that reading one polarity peak value 


is at least a useful guide. The discharge time constant of Cj via the 
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input resistance of amplifier A8 (about 1 second) gives a rough indicat- 


ion of the average signal level. 


sh Siete! THE PDM-MODULATOR TRANSFER CHARACTERISTIC 


The transfer characteristic, from the voltage level at the 
(+) input of the buffer-limiter amplifier A7, to the pulse-duration 
modulation observed at the output of the voltage-comparator A9, is of 


interest from’a system-calibration standpoint. 
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Figure 3.29 The PDM-Modulator Transfer Characteristic 


The detailed circuit schematic, with component values, for 


the preceding sub-sections of Section 3, are combined in Figure 3.30 


on the following page. 
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Sos THE BATTERY POWER SUPPLY 


The detailed schematic of the battery power supply is given as 
Figure 3.31. The following notes cover the design and operating details. 

For availability, weight, and cost considerations, standard- 
size "C"' flashlight cells are used in the battery pack. Twelve cells are 
arranged to provide +9V and -9V relative to ground. Such cells are 
rated at 1200 mA-hours, intermittent use, to a nominal “end-voltage" 
of 1,2 Volts. The total current drain of the E-field sensor and the 
amplifier-filter-modulator system is approximately 26 mA, so that a 
useful battery life of about 45 hours can be expected. Alkaline-cell 
equivalents can of course be substituted for longer service life. 

To retain constant operating voltages for some of the circuits 
(primarily the 1-kHz ramp generator), the battery voltages are reduced 
and stabilized at +6.8V by simple Zener-diode/transistor voltage regul- 
ators. With the low collectoremitter voltages involved, and at the low 
current levels concerned, small-signal transistors (types 2N3904 and 
2N3906) are adequate as series regulators without heat-sink precautions. 
Regulation is sufficient to hold operating voltages at +6.8V +.05V, for 
a battery-voltage run down from 9.5V to 7.4V. 

A battery-test facility is included which makes use of the 
“modulation-level" panel meter. A Zener-diode type 1N5233B suppresses 
6V of the battery voltage, to permit an expanded-scale reading in the 
OVe-eoy LanPeuOreintesest. pliieyy.ok0nm resislor sets meter) cei lection 
at approximately half-scale for the "end-of-life" indication. The +9V 
and the -9V battery packs may be checked individually by means of a 


panel-mounted selector switch. 
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OME THE PULSE-DURATION DEMODULATOR 


During playback of a magnetic-tape recording, the PDM demodul- 


ation involves the following functional circuit blocks: 






Amplifier & 
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Figure 3.32 Demodulator Block Diagram 


The design details and operation of these circuits are covered 


in the following sub-sections. 
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Oise: INPUT AMPLIFIER AND ZERO-CROSSING DETECTOR 


The playback signal derived from the recorded differentiated 
PDM waveform has been sketched in Figure 2.8, but for use by the 
demodulator circuits, this amplitude-variable non-sinusoidal signal 
must be reformed into the best possible replica cf the original PDM 
waveform. A simplified sketch of the circuits involved is shown in 


Figures 3.32). 
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Ficure.3.33 Input Amplifier and Zero-Crossing Detector 


Referring to the above figure, the playback signal from the 
cassette tape recordenms developed across an input-load resistor Ry. 
This 10-ohm resistor provides the load equivalent of an earphone, and in 


the particular recorder used, it was noted that such a load resistor of 
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40 ohms or lower, was necessary to prevent recorder-amplifier oscillat- 
ions. The typical portable recorder can produce a peak undistorted 
stenal of about 35Volts across .uchva load. 

Amplifier Al provides a voltage gain of five, set by the resist- 
ance ratio R3/Ry, so that the +3V peak signal input can slightly over- 
drive the amplifier output (about +14V with the +15V supply voltage). 
Since an output signal of about +10V peak is normally used,full recorder 
output volume is not required, avoiding possible distortion from this 
source. Capacitor C, adds a slight high-frequency emphasis for somewhat 
sharper "zero crossings". 

This precaution against waveform peak clipping in the input 
amplifier is necessary to ensure stability of the regenerated PDM wave- 
form, since time displacements of the zero crossings can occur if the 
waveform is clipped and the duty cycle changes due to modulation. A 
small panel meter (M) is provided so that the peak signal level at the 
output of amplifier Al can be conveniently set to a reference value 
during signal playback. Diode Dj stores the peak (positive) value of the 
signal on capacitor C5, and this voltage is "read out" by the high-input— 
impedance voltage-follower A3. Resistor R7 provides the required meter 
calibration. 

Amplifier A2 is the "zero-crossing" detector, which recreates 
the pulse-duration-modulated waveform, as sketched in Figure 3.34. This 
zero-crossing detection is simply achieved by operating amplifier A2 
ine a feedback resistor,so that its output is always at either the 
(+) ox the (-—) saturation limit. Acting as a voltage comparator, it 


switches state when the (-) input crosses through the nominal zero 
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voltase at the (4) input. In practice, a slight hysteresis is provided 
by resistors Rs and Re to prevent false transitions due to minor noise 
fluctuations of the input signal. The polarity sense of periods Tj and 


Ty have been chosen for convenience in following circuits. 





Output of 
Amplifier Al 








Output of 
Zero-Crossing 
Detector A2 








Figure 3.34 Zero-Crossing~Detector Waveforms 


3.4.2 MEASUREMENT OF PERIODS Ty AND To 


Direct time measurements are made of the pulse-durations Ty 
and T). With a pulse waveform available at the output of the zero- 
crossing detector, with constant-voltage positive and negative limits, 
an operational-amplifier voltage integrator can be used to make a very 
accurate pulse-duration measurement. The typical arrangement is shown 
im simpaitied Lotmeaiuerrcure 3.05(a), Note that dn vthis drawing, the 


integrator-reset and the sample-and-hold switches are sketched as 


80 





81 


simple contact points, although as detailed in Figure 8-50. ssolid-state 


switches (junction field-effect transistors) are actually used. 


Integrator ie Ty | 
ae igs a Sect 
5} Sample & ae ae em ia ets 





(a) Simplified Schematic (b) Waveforms 


EiLeunce 33) Typical Period Measurement 


Referring to the above sketches, note that the negative transit- 
LomneoL! N a7 (the output of the zero-crossing detector) is taken as the 
start of period Tj, and is used to generate an integrator-reset pulse. 
This pulse actuates switch $1 (turns a FET switch ON) which "resets" 
the integrator by discharging capacitor Cj. The negative voltage level 
of Vin (via diode Dj to ensure that only this portion of the waveform is 
related to the measurementof the period Tj) is integrated for the period 
Tj, and develops a positive ouput-voltage Vo- This integrator-output 


voltage is a linear function of time, since Vin is a constant, with the 
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value of V. described by the relation: 


- Ty 

Vo (T}) a (V5 y Vaiode Ri Cy + NG (0) 
which for Ty = 500 wsec, with typical RC component values as noted in 
Hie. 40 ,,and with an ideal reset £0 v¥,(0) = 0, can be evaluated as: 


: = 
V Ci3tds (oxiver) 


ee tose) oe 


Since this V, is directly proportional to Tj, V. becomes a linear 
voltage analog of the period Tj as it varies from 150 wsec to 650 usec. 
The positive portion of Vy, (which represents o>) ds) cheectively 
disconnected from this integrator by the diode D,, leaving the integrator 
input voltage at zero via resistor Ro, so that the output V, remains 
constant. The positive transition of Vj, is used to generate a "sample- 
and-hold" pulse, which momentarily closes switch S9, which charges 
capacitor Cy to the voltage V.. Note that the change in voltage on C9 
may only have to make up for the charge loss during the preceding period 
T, or if may be appreciably different if the value of V, has been changed 
due to the modulation of period Tj. Amplifier A5 is arranged as a unity- 
gain voltage follower with its high-impedance input monitoring the 
voltage on capacitor Cy with minimum charge loss during a period T. 
This sequence of period T, measurement and the storage of its 
voltage analog on capacitor Cy is repeated during each cycle of the PDM 
jest waveform V;,, i.e., at the 1 kHz PDM repetition rate. The output 


of voltage-follower A5 is therefore essentially a "continuous" voltage 
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proportional to the pulse duration Tj. As the duration of Ty changes in 
response to modulation, the output of amplifier A5 "tracks its changing 
value as monitored by the 1000 samples/sec supplied by the integrator A4. 
A small 1-kHz ripple will of course be present on the signal at the out- 
put of amplifier A5 due to a trace of charge loss from capacitor Cy, and 
due to any AT, between succeeding samples. 

Period Ty is measured in a similar manner, but the decoupling 
diode is reversed so that only the naetrente periods ot V7. are used by 
its voltage integrator. 

The total period T is measured separately, and is discussed 
under Section 3.4.4. 

It might be noted that this demodulation method is much more 
effective than the usual sharp-cut-off low-pass filter used for PDM 
demodulation, especially in terms of its transient response. The sample- 
and-hold ripple at 1 kHz is more than 70 dB down from its source 
magnitude as the PDM waveform Vip, and a full-scale modulation shift 
can be tracked in one or two samples, giving a rise time of 1 msec to 


2 msec, with no waveform overshoot. 
Reproduction at the 
Demodulator Output 
of a Recorded 10-Hz 
Square Wave. 


Figure 3.37 Square-Wave Response of the PDM System 


The above chart record shows a reproduced 10-Hz square wave, from modul- 
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ator input, via a magnetic-tape record, and back through the demodulator 
Coma recorder. 
Referring to the detailed schematic Figure 3.36, note that the 
(+T,) and the (-Tj) voltages are summed by amplifier A6, whose output 
voltage therefore represents the signal -(Tj - Ty). Since the measure- 
ment channels for periods Tj and Ty are identical except for polarity, 
Tiptoeevidents toate lt gl =e) representing zero modulation, then the 


output of amplifier A6 will be zero as required. 


324.3 RESET AND SAMPLE PULSES 


All reset and sample pulses are directly or indirectly derived 
from either the leading or the trailing edges of the PDM waveform, as 
regenerated by the zero-crossing detector A2. 

As a specific group of functional circuits, these are detailed 
on a separate schematic Figure 3.38. 

The duration of an integrator-reset or a sample~and-hold pulse 


is nominally 10 usec, which is roughly determined by the time constant 


of the input RC coupling network ( a 150 pF capacitor and a 68K resistor). 


Note that the pulses labelled (P4) and (P5) are derived from the trailing 


edges of the (P3) pulse, rather than directly from the (Pl) waveform. 
This is necessary, since the functions controlled by pulses (P2)Pands( eo) 
must be completed before the functions controlled by pulses (P4) and (PS) 


can begin. All pulses have +14V levels: the positive pulses are used to 


turn ON the type MPF-102 N-channel junction-FET gates, while the negative 


pulses are used to turn ON the type 2N5460 P-channel FET's. 
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3.4.4 TAPE-SPEED COMPENSATION 


The circuits detailed in Figure 3.40 provide the compensation 
for tape~speed variations, using the method discussed in Section 2.4.4 
(the generation of a signal proportional to (T7-T»)/T). The basic wave- 


form relationships are sketched below as Figure 3.39. 
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(a) "Normal" Tape Speed (b) Higher Tape Speed 


Pisures 3.37 Tape-Speed Compensator Waveforms 


(refer tol texetereiubl discussion) 
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The total-period T is measured independently of the T, and To 
measurements. While it is true that T is already known implicitly as the 
sum of T, and Ty> the combined timing errors lead to a value of T which 
is not sufficiently stable to use in tape-speed compensation. 

Referring to Fig. 3.40, amplifier All is arranged as a voltage 
integrator which produces a negative-going ramp of fixed slope. At the 
end of each period T, sample-pulse (P3) gates the value to a storage 
capacitor, with the voltage "read out" as V, by voltage-follower Al2. 
(Fig. 3.39, second row from the top, shows the relevant waveforms and 
their time relationship.) As soon as this sampling is complete, pulse 
(P4) resets the integrator to zero, and the measurement cycle repeats. 

The "current" value of T (the voltage V,) is the input to a 
second integrator Al3, whose output is therefore a positive-going ramp 
whose slope is proportional to the period 1 (the longer the period T, 
the more negative is V,, and the steeper is the output ramp from Al13). 

The ouput ramp from integrator Al3 is applied as a sweep volt- 
age to the (-) input of voltage-comparator Ala the (+) input us held 
at a fixed +5V. The output of Al4 is therefore a waveform with the 
positive-pulse width related to the slope of the sweep voltage at its 
(-) input; the steeper the slope, the less time required for the sweep 
to rise from OV to the +5V reference level, at which time the output 
changes state. If the period T increases, the pulse-width (P7) decreases 
and is proportional to 1/T (3rd and 4th row waveforms of Fig. 3.39). 

The FET switch Q 1 and the integrator A9, constitute a "time- 
amplitude" multiplier. In this case, the amplitude of the input voltage 


represents =(T)-T>) and it is integrated for the time lyjt. Che output 
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of integrator A9, at the end of the (P7) pulse, is therefore proportional 
to (1)-T5)/T and this*value is "sampled" by pulse (P8), stored on a_ hold 
capacitor, and read out by voltage-follower A10. 

Figure 3.39(a) sketches this sequence for a "normal" period (T), 
with the output of integrator A9 being proportional to the area of the 
shaded rectangle (bottom row). Figure 3.39(b) represents the sequence 
during a tape-speed increase, which has slightly reduced the total period 
(T). However, the shaded rectangle now has longer duration but reduced 
amplitude, so that area 1s a constant .(for this particular value of 


inh S To). The design performance is a very close approximation to: 


/ 


yet 2) (y= 1D) 
ak Ts 

This measurement and integration sequence is repeated for each 
cycle of the PDM waveform (Pl), which as already noted, has a repetition 
rate of 1 kHz. The demodulated PDM signal is now a pseudo-continuous 
analog voltage at the output of amplifier A10. This output voltage can 
of course be scaled as desired by selecting the appropriate value of the 
feedback resistor (from the output of Al0 to the (-) input). A simple 
output RC low-pass filter has been used to reduce the trace of 1-kHz 
“sampling” noise, since the recording system used was capable of reprod- 
(Vea Sepa BOE 

The effectiveness of such a tape-speed-variation compensation 
system can be demonstrated by introducing a "synthetic" PDM waveform 
as (Pl) in place of the output of amplifier A2 (Fig. 3.36). This wave- 


form was derived from a standard "function generator", which can provide 
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an asymmetrical pulse waveform, with a swept frequency; the equivalent 
of a regenerated PDM waveform from a tape recorder with a very bad 


tape-speed variation. 


ay Sy ek 


Fi00 msec 


Figure. 3.41 Tape-Speed-Variation Compensation 


For the above test, the pulse ratio was set at Tj = 0.8 T, 
and the frequency sweep was from 1 kHz to 1.2 kHz in 100 msec; a 20% 
frequency shift. Figure 3.41(a) shows the expected change in the demod- 
ulated signal (T,-Tj), (recorded from the output of amplifier A6, Fig.3.40, 
with the D.C. component suppressed). Figure 3.41(b) shows the "compens- 
ated" signal (T1-T2)/T at the output of amplifier A10. The recorder gain 
for (b) was ten times that used for (a), to show the extent of the comp- 
ensation. (The "spike" disturbances in the compensated signal Of Figs 
3.41(b) occur during the frequency-sweep "flyback", which generates an 
equivalent tape acceleration which exceeds the compensator tracking rate. 
This is over 250 times the maximum tape acceleration expected in even 


a poor tape recorder, which might have a 27% speed change at a 10-Hz 


“Flutter rate. 
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It should be noted that this particular characteristic of 
pulse-duration modulation has long been recognized, and under other 
names (such as pulse-ratio modulation) has been used for various 


approaches to tape-speed "wow-and-flutter" compensation [68],[69]. 


sigue POWER SUPPLY 


The demodulator has been designed with an A.C. operated power- 
supply, on the assumption that playback of a signal tape and demodul- 
ation of the PDM signal can be done under “base-station" conditions 
where A.C. power is available. 

The detailed schematic is given as Figure 3.42. This power 
supply provides +15V and -15V relative to a common eround, for use 
by integrated-circuit operational amplifiers. Since the load is 
essentially constant, adequate voltage regulation is provided by 
Zener diodes and series transistor regulators. FET current regulators 
are used to provide constant Zener currents, which improves the 


effective line regulation. 
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325 PHYSICAL-DESIGN DETAILS 


The physical design of this prototype system has been largely 


dictated by the intended portability and component availability. For 


the sake of completeness, the following notes and sketches are included, 


but it should be appreciated that they relate only to this prototype 


system, 


which was intended to evaluate a number of design concepts. Many 


physical-design alternatives are of course possible, and could be 


equally acceptable. 


(a) 


(b) 


The E-FIELD SENSOR assembly is sketched in Figure 3.43. The 
physical dimensions were determined by the convenient size of 
the "Hammond"#1442-15 chassis (10"x6"x1"). This provides a 
sensor-plate area which is large enough to generate a useful 
signal voltage (a larger plate would improve the sional Eo— 


noise ratio) and small enough to fit into an attache case. 


The AMPLIFIER-MODULATOR system was designed to fit into a 
standard black-phenolic "meter case" (6 3/4"x 5 oy ee Oe 
Figure 3.44 shows the panel arrangement used, the printed- 
circuit board location and the battery space. Some care is 
required to mount all of the components of Figure 3.21, 23-50% 
and 3.31 onto a single-sided PC board of the dimensions shown. 
The shield behind the PC board is imperative to prevent local 
pick up of 60-Hz interference by the high-impedance notch- 
filter circuits. Access holes must be drilled through this 


shield to permit frequency and null adjustments. 
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{c) The PDM DEMODULATOR system is also housed in a meter case, with 


the physical layout as shown in Figure 3.45. The circuit compon- 


ents of Figures 3.36,3.38, and 3.40 are mounted on the plug-in 


printed-circuit board. The rectifier and regulator components 


of Figure 3.42 are mounted on a separate PC board secured to a 


sub-panel bracket. The output "banana jacks" are spaced to accept 


standard molded "General-Radio" twin-plug connectors. 


The following are the nominal weights of the system sub-assemblies: 


E-field sensor and 40' of cable 
Amplifier-modulator (with batteries) 


Portable tape recorder (with batteries) 


Attache case 


Field-portable system 


PDM demodulator 


R 
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3.6 THE MAGNETIC-TAPE RECORDER 


The portable magnetic-tape recorder used with this system is a 
"Realistic'' Model CTR-20 (listed as "Radio Shack" Cat. No. 14-873). The 
nominal frequency response is given as 60-10,000 Hz, and has been 
adequate for the required time resolution of the PDM pulse edges. 

During recording, the output signal from the amplifier-modulator 
is introduced (via a shielded patch cord) into the "AUX" input of the 
tape recorder. The optimum recording level for the differentiated PDM 
signal appears to be at approximately -10 dB. 

During playback, the recorder output signal is taken from the 
"FAR" jack, and introduced (again via a shielded patch cord) into the 
input jack of the PDM demodulator. Proper playback level (to prevent 
input~-amplifier overload) is determined by observing the signal-level 
indicator on the demodulator panel. 

As with all magnetic-tape systems, some problems related to 
signal "drop-outs" have been encountered. If proper zero crossings are 
not developed in the demodulator, valid Ty and T, pulse widths cannot be 
established, and the output signal will be invad d ebOr one Ory more 
samples. These effects are illustrated in Figure 3.46. The recording of 
Figure 3.46(a) originated from the first few feetror tape of a well- 
used cassette. A visual check showed small "Norease'" marks and "burnished" 
spots, which were presumably made by the drive-spindle rotation on a 
momentarily-stalled tape, probably during repeated replays and rewinds 
of this section. Figure 3.46(b) shows the signal derived from a later 


section of the tape, but with a somewhat low setting of the recorder 


volume control, leading to marginal control of the zero-crossing detector. 
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(c) 


(a) Poor section of tape (worn spots). 
(b) Good tape, but low playback volume. 


(c) Good tape and proper playback volume. 


Figure 3.46 Signal "Drop-Out'"' During Tape Playback 


The recording of Fig. 3.46(c) was from the same section of tape as was 
Piews.460(b),, bUL with the volume-control setting increased until 
proper switching of the zero-crossing detector was observed. 

Any good quality magnetic-tape cassette appears to be adequate, 
providing due care is exercised to prevent the physical blemishes noted 
above. 

It should be noted that the presence of a few drop-outs in an 
otherwise good recording need not destroy the value of the record. The 
drop-outs rarely last for more than a few samples and typically produce 


output-signal discontinuities which are recognizable; the tape-signal 
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OCR. 


drop-out results in an immediate period-measurement error (oleae alo.Or eh 
which is sampled and produces an output-signal discontinuity which has a 
rise time limited only by the output RC filter and the recorder writing 
speed. The observed drop-out artifact is typically of several-volts 
magnitude and 2-4 msec rise time, which with a chart speed greater than 
a few inches/sec has a visually-different characteristic from normal 

ELF signals; thay can therefore be ignored, or the playback can be 
rechecked with different adjustments. If the signals are digitized for 
computer data reduction, such invalid drop-out spikes can be detected by 
a simple algorithm related to their excessive dv/dt, and the invalid 
points can be dropped and replaced in some acceptable manner (linear 


interpolation, a running average, GECKO NR 
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3.7 SYSTEM TRANSFER CHARACTERISTICS 


The transfer characteristics of the system sub-assemblies have 
been noted in their respective sections. The following notes cover the 
over-all system characteristics, from incident E-field to the final 
demodulator output record. 

The system frequency response is given in Figure S41. Lie) OUahzZ 
notch width is less than 1 Hz, and the dotted response above 200 Hz is 
an indication that the effective amplitude response becomes datiien lito 


define as the frequency approaches the 1 kHz PDM repetition-rate. 


0,1 1 Baie 1000 4, 


dB 








Figure 3.47 System Frequency Response 


Although such a frequency-response curve implies a sinusoidal 
response, this is not true as waveform degradation occurs in a sampled- 
data system. Above 150 Hz, waveform reconstruction is dependent on the 


available number of samples per cycle provided by the PDM “carrier” 
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system. Figure 3.48 shows this loss of waveform reproduction as the 
number of samples per cycle drops from 10 samples/cycle at 100 Hz, to 
only 2 samples/cycle at 500 Hz. The output RC filter does provide a 
slight smoothing of the sampling steps without excessive slowing of the 


pulse rise time. 


Chart Speed 


Freq. (Hz) _opeed 
100 VAVAVAVA wevaval 10 msec/cm 
200 LUAU AWA Ey AV 5 msec/cm 


00 [WAN VV 
400 AVA V AAU NACA NAV UAV TAU " 
500 DD ARR DLP LI LDL LID LILIIING " 


Figure 3,48 Sampling Rates and Waveform Reconstruction 


The system pulse response is of course degraded by the low-freq- 
uency roll-off below 5 Hz, which in the time domain implies a time cons- 
tant of about 30 milliseconds. Figure 3.49 shows the system pulse resp- 


onse from incident E-field signal, through the PDM-modulation process, 
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magnetic-tape recording and playback, 


ing. 


(a) | os i. . 


(b) 3 Pe ars Sag, ae 


(d) 


(a) Input waveform: 10 msec pulses. 

(b) Demodulator output: vert.= 5V/cm, hor.= 20 msec/cm. 
(Input E-field pulses = 10 mV/m; Atten. = 0 dB) 

(c) Demodulator output (scale factors as in (b) above) 
(Input E-field pulses = 100 mV/m; Atten. = —20 dB) 

(d) Demodulator output (scale factors as in (b) above) 


(Input sensor shielded; Atten. = -20 dB) 


Figure 3.49 System Pulse Response 


to demodulation and chart record- 
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Figure 3.49(a) is the voltage waveform which establishes the 
incident E-field at the sensor plate. Figure 3.49(b) displays the full 
system sensitivity, with a 10 mV/m E-field pulse of 10 milliseconds 
duration recorded with 0 dB attenuation. Figure 3.49(c) reproduces a 
similar pulse but with an amplitude of 100 mV/m recorded at -—20 dB 
attenuation. Figure 3.49(d) shows the system noise level, at an attenua— 
tor setting of -20 dB but with a zero input signal (an electrostatically- 
shielded sensor assembly). Note that with a positive voltage pulse 
applied to a calibrating plate mounted above the sensor plate, the 
incident E-field is negative. The system output has been arranged so 
that the strip-chart recording polarity agrees with the polarity Of Ehe 
incident E-field: 

Figure 3.49 can also be used to calculate an over-all transfer 


characteristic, as shown in Figure 3.50. 
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Figure 3.50 Input-Output Transfer Characteristics 
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Figure 3.50 provides an over-all system calibration: the demod- 
ulator output voltage versus the amplitude of the incident E-field signal 
(assuming the frequency is within the 5 Hz-200 Hz pass-band) for various 


settings of the amplifier-modulator attenuator switch. 
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CHAPTER 4 


CONCLUSION 


4,1 TYPICAL FIELD RECORDS OF ELF SIGNALS 


As noted in the introduction to Chapter 3, the availability of 
a Hewlett-Packard Model 5480A Signal Analyzer has permitted the reprod- 
uction of field records with a wide choice of vertical sensitivity and 
horizontal time scaling. The resolution however is limited by the 1000 
samples/scan, and in some of the following illustrations of expanded- 
scale X-Y recordings this may imply only 50 samples/cm. As a result, the 
apparent fine structure observable is in some cases simply quantization 
noise and must be interpreted with some caution. 

The following pages show typical records of natural ELF signals. 
The vertical scale is given in terms of mV/m E-field amplitude, and 
the relevant horizontal time scaling is noted. A few features OfmInveres. 
have been marked, and selected eerione have been repeated with an 
expanded time scale. 

These particular records come from a series made at approximately 
0800 hrs M.S.T. May 17, 1975, in Mayfair Park, Edmonton, Alberta. This is 4 
river-valley municipal park, located close to the centre of the city, 
with the actual recording site being an open grassy area some 100 yards 
away from the nearest roadway and from lighting standards. The sky was 
cloudless and only a light breeze was blowing. No data are available 


regarding continental lightning activity for this period. 
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(a) Vert: 40 mV/m/cm Hor: 1 sec/cm 
(b)(c)(d) Vert: 30 mV/m/cm Hor: 1 sec/cm 


Scaling Data: 


(e) Vert: 12 mV/m/cm Hor: 0.5 sec/cm 


Features: £1,f5,f3,f4.---typical bursts of 1 or 2 cycles, with a 
period of about 150 msec., initial deflection may be 


either positive or negative. 


nA ee burst of oscillations at about 13 Hz, presumably 


related to the n=2 mode of the Schumann resonances. 


Figure 4.1 Typical ELF Field Recordings (Sample #1) 
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Bo ee 
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Expanded as (c) 
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Scaling Data: (a) (b) Vert: 30 mV/m/cm Hor: 1 sec/cm 
(c) (d) Vert: 12 mV/m/cm Hor: 0.5 Sec/cm 


Features: f,....short bupetasimilarscosty,.7.0leeie. Ui Wee Pel ckae, 
f>...."sawtooth" pattern at about 4 Hz, also noted in 
other portions of this record. 
f3....short burst Sinulan tomb, buteat about 13 HzZe 


f,,f5....bursts of several cycles at about 5 Hz. 


Figure 4.2 Typical ELF Field Recordings (Sample #2) 





The preceding ELF recordings (Figures 4.1 and 4.2) have of course 
been selected because of their "features" of interest. Minute-to-minute 
changes in the character of the ELF signals are evident, primarily as 
changes from relatively-quiet periods (such as in Figures 4.1(a) and (b)) 
to periods displaying various types of noise-burst activity (such as in 
Figure 4.2(b)). 

The noise bursts recorded in Figures 4.2(c) and 4.2(d) display 
several characteristics which might be noted. Features such as (£4) are 
frequently observed; typically little more than one cycle with a period 
of about 150 msec (26.7 Hz) and presumably related to the n=l mode of 
the Schumann resonances. Feature (£3) is a similar burst but at about 
13 Hz (presumably the n=2 mode). In each case, the resonant-mode oscill- 
ation is abruptly terminated by an "over-riding" signal having a distinc- 
ly different characteristic; in the case of (£3) the new signal (f£,) 
consists of several cycles of a damped 5 Hz oscillation. Feature (£5) is 
a rough "saw-tooth'"' waveform with a fundamental frequency of about 4 Hz, 
and this tendency towards a saw-tooth pattern is noticeable in other 
parts of the same recording. 

The recording samples on the following pages (Figures 4.3, 4.4 
and 4.5) were made at the same Mayfair-Park site as those of Figures 4.1 
and 4.2, but several months later; at 1300 hrs M.S.T. September 26,1975. 
Some of the noise bursts on these later records have a greater complexity 
and generally a somewhat greater amplitude than those of the earlier 
recordings. There has been insufficient monitoring to offer any comments 


as to whether these are diurnal or seasonal differences. 
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Scaling Data: (a) Vert: 30 mV/m/cm Hor: 1 sec/cm 
(b) Cc) Vert: 12 mV/m/cm Hor: 0.25 sec/cm 


Features: iq. . .0sciliations at approximately 2/7 Hz, superimposed 
on a low-frequency (5 Hz) disturbance. 


fio. -weeeveral oscillations at. about 22 Hz. 


f£3....oscillations at about 22 Hz, superimposed on other 


disturbances. 


Figure 4.3 Typical ELF Field Recordings (Sample #3) 
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Scaling Data: (a) Vert: 30 mV/m/em Hor: 1 sec/cm 
(by Gc) Vert: 12 mV/m/cem Hor: 0.25 sec/em 


Features: £,..--oscillations atabout uo GHz. 
fo....disturbances with components at 5.7 Hz, 
Jeg suze Aowie and 2020z), 
f,....disturbances at about 13 Hz, but less coherent 


than the burst noted as fi. 


Figure 4.4 ‘Pilypleal ELLY Field Recordings (Sample #4) 
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Scaling Data: (a) (@)(c) Vert: 12 mV/m/em Hor: 0.25 sec/cm 


Features: f,.---damped oscillations at about 21) Hz. 


f5....oscillations at 15 Hz superimposed on other 


disturbances. 


f....-negative impulses at a 6.7 Hz repetition cate. 


f,-...oscillations at about oz HZ: 


Figure 4.5 Typical ELF Field Recordings (Sample #5) 
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In Figure 4.3(b), feature (f,) appears to be a lightly-damped 
burst of oscillations at 27 Hz (the n = 4 mode of the Schumann reson- 
ances) superimposed on a few cycles of a 5 Hz disturbance. A common 
excitation source might be presumed because of their time juxtaposition. 
In Figure 4.3(c), feature (fo) Pcons tists of damped oscillations at about 
22 Hz (perhaps the n = 3 mode), and feature (£3) is a longer burst at 
22 Hz superimposed on a low-frequency disturbance, andspurldineg upsto a 
higher amplitude for a few cycles. 

In Figure 4.4(b) feature (f;) consists of a burst of damped 
oscillations at about 13 Hz (the n = 2 mode). Similar components can be 
seen in feature (£5) of the same recording but incoherently mixed with 
disturbances at several other frequencies. Further traces of the 13 Hz 
oscillation are evident in recording-strip (c) feature eae 

In Figure 4.5(a), feature (£4) is a damped oscillation at about 
21 Hz (the ni= 8 mode). In recording-strip (b), feature (£5) consists 
of 15 Hz oscillations apparently excited by the larger negative impulses 
which occur at a repetition rate of about 6.7 Hz. These impulses do not 
appear to be fundamental-mode oscillations since they are not sinusoidal, 
but the uniformity of spacing and the similarity of their amplitude and 
shape would imply a common origin. Recording-strip (c) includes bursts 


of oscillations (such as feature Cee. at about 32 Hz (the n = 5 mode). 
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4.2 COMMENTS 


Several comments should be added to the material covered in the 
preceding chapters, and these comments relate to two specific themes: 
circuit design details which have been tolerated in the prototype system 
but are quite evidently not optimum, and secondly the problem of system 
transient responses and possible signal artifacts. 

From the standpoint of further circuit development, the relevant 


points might be itemized as follows: 


(a) The newer high-grade low-noise FET operational amplifiers (such 
as the Burr-Brown 3523 series, or the Teledyne-Philbrick 1421 
series) should be evaluated for use as the electrometer/cable- 


drive amplifier in the E-field sensor assembly. 


(b) The battery drain of the amplifier-modulator module could be 
considerably reduced by the use of low-current-drain operational 
amplifiers (such as the RCA CA-6078A) in several locations. 

Such units could be used in the linear-amplifier stages although 
their slew-rate limitations would probably rule out their use in 


the pulse circuits. 


(c) The temperature stability of the prototype system has not been 
evaluated. Voltage-offset drifts are a minor concern since the 
closed-loop gains are relatively low, and integrator drifts lead 
only to D.C. shifts which are zeroed out during playback and 
the generation of a chart recording. The system has been used in 
the field over a temperature range of O° 6-125 "Co with no) evidence 


of trouble except an occasional sign of null shift in the 60-Hz 
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notch filter. Metallized-film polycarbonate capacitors with a 
temperature coefficient of less than: 100ppm/°C were used in 
obviously critical circuits, but a complete review of the 


system temperature stability would be useful. 


(d) It would be interesting to explore the added capabilities of 
a "stereo" portable recorder. A second signal with a pass band 
from 200 Hz - 10 kHz could be extracted from the electrometer 
amplifier, and this could be directly recorded as an analog 
signal on the second recorder channel. On playback, time- 
correlated ELF and VLF signals would then be available for 


recording on a 2-pen system. 


With regard to system transients and potential signal artifacts, 
the general problem has already been discussed in Section 273. tOPeLuer 
with examples of pass-band transients generated by filter step and 
impulse responses. However, a comment should be added in regard to the 
apparent lack of concern by investigators in the ELF field regarding 


such instrumentation problems. 


In reviewing the rather extensive literature on ELF investigat— 
ions, it was noted that there was no documentation of instrumentation 
step or impulse responses, no mention of the type or order of the low and 
high-pass filters used, only an occasional sketch of the frequency 
response of a system and never a full Bode plot of magnitude and phase 
versus frequency. Only a single mention was found of the possible degrad- 


ation of waveforms by notch filters (/O)'= Unisslackeok concern regarding 


system transient response is all the more surprising in view of the 
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impulsive nature of typical ELF signals, and the fact that some aspects 
of theoretical work have been based on the analysis of waveshapes of 
recorded transients. 


Figure 4.6 is an illustration of this point. 


250 usec 
a 
ee Input Pulse 
Output From a 
(b) 


Low-Pass Filter 


(c) (eregers With an Added 


k—>+ msec 60-Hz Notch Filter 


A "Slow-Tail" 
ELF Waveform 


<5 msec (after Hughes [71]) 


Figure 4.6 Filter Artifacts and a "Slow-Tail" Waveform 


In the above Figure, waveform (a) is a 250 usec input pulse, 
waveform (b) is the resultant output from a 2nd-order Butterworth low- 


pass filter which has a 500-Hz upper corner-frequency, and waveform (c) 
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shows the added distortion if a passive Twin-T 60-Hz motch filter is 
included in the signal path. Wavetorm (4) is a sketch of a typical 
ELF "slow-tail' waveform used in the study of far-distant lightning 
discharges [71],[72]. It is evident that one needs only select the 
proper order low-pass filter with an appropriate corner-frequency, 
and the filter impulse response would match the "Slow-tail" waveform. 
A somewhat related concern pertains to published recordings 


such as those reproduced in Figure 4.7. 
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TYPICAL ELF ATMOSPHERIC NOISE 
Wide-band and narrow-band ELF noise from Florida (Feb. 1968). 











Figure 4.7 Typical ELF Atmospheric Noise (Evans and Gritrithes (7 oi) 
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The implication is that the wideband noise contains the narrow- 
band components illustrated. However, the impulsive nature of the wideband 
record is evident, and a transient pulse a few milliseconds in duration 
obviously does not have a 1/2-sec train of 20-Hz oscillations as a comp- 
onent. The further implication that low-frequency components maintain 
phase coherency over dozens of cycles is also unfortunate, and does not 
agree with reality; waveform coherency for more than 4 or 5 cycles is 
rarely observed. The narrow-band recordings of Figure 4.7 simply portray 
the impulse responses of the high-Q sharply-tuned filters. Any attempt 
to use such equipment in making estimates Of spectral densities would 
seem to be rather futile. 

An additional point in this regard is that published spectral- 
density estimates [15],[17],[22] have probably been based on the classic— 
al method which compensates for the instrument frequency response after 
the time-domain data has been transformed to the frequency domain, since 
multiplication is then involved instead of convolution. Unfortunately, 
there is no evidence that anything except the steady-state frequency 
response of the system is ever used for this correction. From a circuit— 
theory standpoint, it has been pointed out [74] that a short-burst of 
oscillations is equivalent to multiplying a steady-state sinusoid by a 
@ating pulse. in the frequency domain, this requires the convolution of 
the instrument's steady-state response with some sinc function, to 
derive its "impulse" frequency response. Since this transient frequency 
response may be considerably different from the steady-state response, 


the published spectral-density estimates may be considerably in error. 
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The Extremely-Low-Frequency (ELF) spectrum of electromagnetic 
waves (3 Hz—3 kHz) has a number of special characteristics. The back- 
ground noise in this range is generated primarily by global lightning 
aetivity, and tne resultant disturbances may propagate to great distan- 
ces in the Earth-ionosphere cavity (via a TEM guided-wave mode) since 
the attenuation at such low frequencies may be less than 1 dB/1000 Kn. 
Cavity resonances are identifiable at AbouUteouHa, 14 HZ, 20 HZ, <5 5 
the so-called Schumann resonances. The extensive earth penetration at 
ELF wavelengths has been used for geophysical prospecting, and penetrat- 
ion is sufficient even in sea water to permit signal transmissions to 
submerged submarines. 

Measuring and recording systems for ELF signals are typically 
semi-permanent large-scale facilities. Directional H-field sensors are 
heavy, bulky and must be specially installed to be free of vibration in 
the Earth's natural magnetic field. Vertical E-field components are 
usually monitored by omni-directional elevated-whip or ball antennae. 
Extensive back-up instrumentation is typically required, with signals 
being recorded on computer-grade magnetic-tape systems for subsequent 
analysis. 

This thesis is concerned with the system design of a portable 
receiver-recorder for special field investigations. The nominal pass - 
band is 5 Hz — 200 Hz to cover all Schumann-resonance phenomena with 
acceptable waveform reproduction. Battery operation and portability 


has been a prime consideration (13 lbs total weight) and a cassette 
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magnetic-tape recorder provides the data-storage medium. 


The system design has been given in detatl. with all relevant 


design equations and circuit schematics for: 


1) A "charge-induction" type of E-field sensor, with an electrometer 


input and a cable-drive amplifier. 


2) A semi-remote (40 ft from the E-field sensor) variable-gain 


3) 


amplifier system. The frequency response is shaped to minimize 
pulse-waveform distortion, and a modified 60-Hz notch filter 
has been included which does not generate pulse transients. A 
Pulse-Duration-Modulation (PDM) system is included to permit 


the recording of ELF signals on magnetic tape. 


A PDM demodulator system permits recovery of the ELF signals 
from the magnetic tape. Special circuits are included to compen- 
sate for tape-speed variations during playback. Output of the 
regenerated ELF signal is at a suitable voltage level for chart 


recording or analog/digital conversion. 


Reproductions are included of typical field-recorded ELF 


signals, and these records ane at Least’ equal ain quality to those in 


the published literature on ELF research. It is noted that signal 


artifacts are especially likely in ELF instrumentation, due to the 


impulsive nature of ELF signals and the usual transient responses of 


signal) filers. 


The detailed derivation of the modified 60-Hz notch filter is 


given in Appendix A, and illustrates how transient and steady-state 


responses can be separated to ensure a good pulse response. 
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APPENDIX A 


THE MODIFIED TWIN-T 60-HZ NOTCH FILTER 


A.1 THE PASSIVE TWIN-T NOTCH FILTER 


A typical twin-T network, with symmetrical elements and a 
Wee even loadecrecnowne in Pieureral (a), and its frequency response is 


sketched in Figure A.l(b). 








(a) (b) 


Figure A.1 The Passive RC Twin-T Notch Filter 


The voltage transfer function Lote icine Al (a) can be written’, 


in standard Laplace notation, as: 





V (s) = H (s) = pee es a tu Gees = 
Vie) o s2 + 2(2 + R/Ry)wos + (2R/Rz, + 1)us 


where w, = 1/(RC). With negligible loading on the filter output, €-2-,5 
4f the notch filter feeds into a high-impedance voltage Followers — hyo > 


and-the transfer function reverts leon 


32 ae we 


PAG) 2 Goeeenie SESE (A-1) 


s2 + 4uwos + w4 
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Pie Oot such) awitter is. low (.25) and for a 60-Hz notch frequency 


this leads to a 3-dB-down notch width of 240 Hz. 


The poles of this filter are real, and the pole-zero locations 


are sketched in Figure A.2. 
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Figure A.2. Poles and Zeros of the Passive Twin-I Filter 


The impulse response of the passive twin-T notch filter can be 


derived from Eqn. (A-1): 


D 2 
Ee Pee Saige oe 
Ho (s) a2 ies + ws 
Fo lw 4. 3w 
> + oo SRO te ES eel 
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and returning to the time domain via the inverse Laplace transforms, 


gives the impulse response: 


Heche [eos Bie eee obese Petal uta) 


which for £. = 60 Hz, t> 9,, and t = milliseconds, gives: 
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The form of this impulse response is sketched in Figure iN 
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Figure A.3 Impulse Response of the Passive Twin-T Filter 


The pulse response can be evaluated by using a unit-amplitude 
input pulse of width (1): 
h(t) 
c(t a= at a eer) 1 
giving, eat (1-e7T$) | - 


S age ie 


The ouput response of the passive twin-T filter with this input pulse, 
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and reverting to the time domain, with £,=60 Hz, and t=msec, the pulse 


response is given by: 


vat) ell ae WeloCemes te sees SS) Jute) 


Si reeieces oD - aa Ont (ET) (een) 


which produces the following waveforms for typical values of the pulse 


width (t): 
ae ae |--—_1—_+ 
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Figure A.4 Pulse Response of the Passive Twin-T Filter 
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Are THE HIGH-Q TWIN-T NOTCH FILTER 


To increase the effective Q of a twin-T notch iulLter.. an active-~ 
network version can be arranged as shown in Figure A. 52), with its 


typical frequency response shown in) Figure A.5(b). 





(a) (b) 


Figure A.5 High-Q Active Twin-T Notch filter 


In this case, if the voltage-follower gain is very close to unity, some 
deliberate network loading by resistor Ry, is usually necessary to reduce 
the gain and maintain stability at the high 0 value desired. 

The voltage transfer funetion for the circuit of Fig. NE EW 6 


assuming a voltage-follower gain Of Unity, will be. 


V2 (s) s2 + wé 


and for the pesietance bya the poles and zeros coincide on the jw 
axis. 
For typical values of components for £,=60 Hz (On=2 tio aes > 


let R=80Kohms, C=.033uF in Fig. A.5(a) above, the gain is readily cont- 
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rolled by varying the resistance of Ry, which can therefore be used to 


control the Q and hence the notch width. 





Figure A.6 Root Locus for the High-Q Twin-T Notch PabLcen 


The locus of the positive complex root 4s shown in Figure A.6, as the 


gain is varied in terms Ofethesload fesistor Ry; the poles become comp lex 


for Ry,>33Kohms. From Eqn. (A-3) the Q can now be expressed in terms of 


this loading resistor Ry; as: 


q =th= Ry 
2R 160Kohms 





aoe, awe eye = 500Kohms (and R = 80Kohms as already noted), the filter 
Q = 500/160 = 3.1, so thactfor to = 00 Hz the notch width has been 
reduced to 60/3.1 = Salo this, notch characteristic is sketched in 


Figure A.5(b); the presence of the pole is evident on the high-frequency 


ec Cee ko. 
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The impulse response will now reflect the presence of the 


complex poles. From Eqn. (A-3), 


s2 + w2 
Ap llye. Saye ree REE le eee 
1 s2 + 2(R/Rz,)wos + (2R/Ry + 1)wo 


whiten Lor Ry = 470Kohms (the closest standard resistor value to the 


500Kohms), gives in terms of partial fractions: 





aes Gli (286) Goon) | 


= = 2 eS SS 
BIRge, ORM Me z F ci7ue)2 + (s15ue)* (+ -17u.)? + (1.1500) 


anduror ft. = 00-82, = = milliseconds, and t>0, the time-domain impulse 


response becomes: 


pie eomeanes 94 = (132 cos 4 otat 114sin.43t) (A-4) 
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The Impulse Response of a High-O Twin 
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The pulse response can be derived by using the same approach as 
for the passive twin-T filter, i.e., the input pulse is assumed to be of 


unit amplitude with a duration T, giving: 


Y,(s) = X(s)H, (s) 
7 Geni) ; g2 + we 
7 Ss 6? + 2(R/Rz)¥os + (2R/Ry, + 1)w2 


and again with resistor ae 470Kohms, the output in terms Offa, partial 


fraction expansion becomes: 


s mie .258(s +.170.) 2264 (1. U5.) 
ee aH tS | 24 2s 9 = © 
tee | =F Gedciyug)? + (1.18u)2 | Gian)? + GQ-150.)° 
Returning to the time domain, with f. = 60 Hz and t = milliseconds, gives 


the pulse response as: 


yz (t) = .74 Pema (. 258cosmaste= 264sin. 43t)| w(t) 


a |.74 Pro elt) (258eos.43(t—t) — 264sin.43(t-1)]] w(e-) 
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which produces the following waveforms for typical values of pulse width: 
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Figure A.8 The Pulse Response of a High-Q Twin-T Filter 
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A.3 THE 60-HZ PEAKED AMPLIFIER 


As moted in Sections A.1 and A.2, both the low-Q and the high-Q 
versions of the twin-T notch filter display poor transient responses, 
although the steady-state frequency response can be shaped as required. 
One approach to an improvement of the pulse response is to use the 
high-Q twin-T network to form a 60-Hz peaked amplifier, and use its 


steady-state output to "Buck out" the interfering 60 Hz in an auxiliary 


amplifier. 
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Figure A.9 The 69-Hz Peaked Amplifier 


For the peaked-amplifier section of Figure A.9, using the high-Q 


transfer function H;(s) of Eqn. (A-3), 


15, (6) ee ee Se) ee ees eee 
Ie Ry | 1 + By (s)Rp (Ry + R5)/ (R3R5) 
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Using typical resistance values (as enoted in Fig. 3.21), a 47Kohms , 


Ry = 47Kohms, R3 = 10Kohms, Ry = 33Kohms, and Rs = 2.2Kohms , 


V5 (s) 7 lt 
= if, E 


2. y 
Ge sr (OF 
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<4 
fe 
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VY 
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2 2 
aus Sei eee Sat he oA). oS 
Hy (s) Se E70 Sues Olus (A-6) 


Evaluating |H, Gu) | gives the frequency response: 





Figure A.10 The Frequency Response of the 60-Hz Amplifier 


The pole~zero locations forthis peaked amplifier are of interest. 


From the transfer function Hy (s), (Eqn. A-6), it can be shown that: 
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It 
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The zeros are at: Ss cla thsy iis 


HT} 


The poles are at: Ss (2700252 71.005). 


The complex pole and zero in the upper-half S plane, with an expanded 


scale, are shown in Figure Nol 
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Figure A.11 Pole-dZero Locations for the 60-Hz Peaked Amplifier 


Note that the pole is safely in the left-half plane for stability, but 
its nearness to j = w. indicates the very high effective 0 (@=2220 8). 
The zero location gives a slightly higher attenuation on the nven—£ red] 
uency side of W., as is evident in the frequency—response curve of 


Figure A.10. 


A.&4 THE MODIFIED 60-HZ NOTCH FILTER 


As is well known, if H(s) is a band-pass function, then [1 - H(s)] 
isthe equivalent band feqection. Ihils (is the effective response of Fig- 


ure A.12 if certain precautions are taken: 
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Figure A.12 Block Diagram of the Modified Notch Filter 
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The voltage transfer function of Figure A.12 is: 


tee) Ds y af gle KH (s) 
Vans) Eee mR Se SakeT | 


Substituting for Hj(s) from Eqn. Kegeeandelettimeeresistor “Rj7= Bo; 





gives: 
2 Ve 2 
H3(s) = - [1 - mole ya ee 
3 \e2 +.0045u,s + 1-013 
R 
Mertings hyo iho =2(.013), rearranging, and ignoring the sign inversion 
5 


of the summing amplifier, the transfer function can be rewritten as: 


2 20045 ae] ee) 2 
S ( ‘pete a ceat te Ky WS 


Pe on sunce al. Olu 


H,(s) = |H3(s)| = C1-Ky) 


By taking the precaution that resistor R3<Rq, (and) that tne nominal gain 


of Ho(s) was unity at £,.), a point can be found on potentiometer P, such 


R 
that its fractional setting = Ko ae For this setting, the 2nd term 
2 


in the bracketed numerator can be made equal to zero, leaving: 


2 1.005w | 

S + ° 

i band . SOR WiC tN a ae es ae _— 
Hy (s) 2! \ SEN, Tee uallg ie oun (A-7) 


In terms of the pole-zero locations: 
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Figure A.13 Pole-Zero Locations of the Modified Notch Filter 
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Note in Figure A.13 that the complex pole location remains the same as 
in Figure A.11l but the zero has been restored to the jw axis, very close 
to the desired j = W.- In practice of course, the frequency-trim adjust- 
- ment in the twin-T network permits a precise location of the frequency 
TUL 

Figure A.14 shows the character of the high-Q notch, and since 
the Q is now over 200, the bandwidth has been reduced to about 60/200 = 
0.3 Hz. (Lf such a notch width is censidered too marrow, the Q’can be 
reduced as required by reducing the resistance valleron Ry in ig: ING) 


which manipulates the pole location asin Figure AL Ot) 





Figure A.14 Frequency Response of the Modified Notch Filter 


While the frequency response of the modified notch filter is 
simply that of a high-Q filter, the impulse and transient responses 
are much different. The impulse response can be derived from Eqn. A-/ 
as follows: 


52 + 1.0053 
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Soi) (oa. Lae 01wé 
: Ss. + .0045wos + 1.01w3 


Hy, (s) 


on a ee ee es eee) 005) ae 
""\(s +.002w.)* + (AN005t6) = (s +.002w.)% + (1.005w.) 


which gives a time domain response, wiehetoes OU) HZ, © = msec, evaisl (8 ea 1e}n, 


hy, (t) = —- 1.67 e0007t (cos.38t + 2 0 s5in.38t) (A-8) 


4, (t) 





Figure A.15 The Impulse Response of the Modified Notch Filter 


The impulse response now reflects the changed pole location; its proxim- 
ity to, the J’ AxiSe produces Very little damping of the sinusoid, but since 
the pole is very close to the zero the amplitude of the sinusoid is only 
about one tenth that of the straight high-Q twin-T filter (as shown in 
alters UNAT OM 

The pulse response of the modified notch filter, with aeunit= 


amplitude input pulse of duration T, Can now be derived from the transfer 
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Y4(s) = X(s)H, (s) 
Teens S2) PeeO0she 
ee 09) a 
iss 4 sb 0045u6s + L.0line 
pes 799 .O1ls -— .0045w 
= SETS |) Se 2 OWES) Soe 
Coke »( eC enas o ai 
| Gere) 299") N1i(s +.002200) 044 (1.0050) 
s (2900220) (00500) Sats +.0022W0)* + (1.005w.) 
Returning again to the fame domain withers =)0U) Hz and ty = imsee,) cne 


pulse response becomes: 


yy, (t) = 99 Hear 0008E ( Oicos.38t - .004sin.38t) u(t) 


- 99 + ee ove (tar) (20lcos.23(t-1) 7 J004sin.38 (tr) u(t-T) 


which is illustrated by the following test results, for several input— 


pulse widths: 
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Figure A.16 The Pulse Response of the Modified Notch Filter 
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It is apparent from Figure A.16 that the pulse response of the 
modified notch filter is a considerable improvement over that of the 
standard high-Q twin-T filter (as shown in Fig. MODs 

From Eqn. A-9 for the time-domain response y,{t)> it can be 
seen that the (.02) cosine coefficient represents a nominal 2% 60-Hz 
ripple on the pulse response, and its exponential damping is negligible. 


This ripple originates as the "error" signal at the output of the 60-Hz 


peaked amplifier after each excitation by a pulse edge; it is not related 


to the 60-Hz content of the incoming signal. (If desired, the level of 
this ripple could be reduced by a higher loop gain and further optimiz- 


ation of the 60-Hz peaked amplifier.) 


A.5 PHASE RESPONSES 


It is of some interest to compare the phase responses Cite 
various versions of the notch filter. From Eqn. A-1l for the tCranusber 


function of the passive twin-T network: 


0 Eee 


iB Coe 


HoGe) = T=GJa.)? + Guo.) 


The phase angle is therefore: 


1 = (w/o) 


6 (a cao ir (A-10) 


In a similar manner, from the transfer function Eqn. A-3 for the high-0O 


twin-T notch filter, the phase response can be expressed as: 
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and from Eqn. A-/ for the modified notch filter transfer function, 


— a1 |_20045w/ we (A-12) 
acs en ~ cord 


The phase responses given by equations A-10, A-11 and A-12 are plotted 


below: 
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Figure A.17 The Phase Responses of the Various Notch Filters 
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A.6 EFFECTIVENESS AND LIMITATIONS 


The effectiveness of the modified notch filter is shown in 
Figure A.18. The input signal consisted of a 13-Hz square wave plus a 
superimposed 60-Hz sine wave, of equal peak-to-peak amplitudes. The 
recorded output is shown (a) with the notch filter switched OUT, and (b) 
with the notch filter switched IN. (The square-wave droop is due to the 


capacitive coupling into the amplifier section which precedes the notch 


Sir eer.) 
(a) yn Notch Filter OUT 
(b) Notch Filter IN 


Figure A.18 The Effectiveness of the Modified Notcisieii ter 


The limitations of this modified notch filter are related to the 
very high "internal Q" of the 60-Hz peaked amplifier. The approach to 
steady-state operation of the system is determined by the exponential 
term of Eqn. A-9, which in units of seconds is exp(-0.8t). The time 
constant is therefore 1.25 seconds, and to settle to within 14 of the 


steady-state value will require about 5.8 seconde. this is illustrated 
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in Figure A.19, which was generated by abruptly eating a constant-ampli- 


tude 60-Hz signal into the notch-filter system. 
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Figure A.19 Settling Time of the Modified Notch Filter 


Under the usual field conditions, the 60-Hz interference is 
essentially constant, so that after the initial turn-on settling time 
Ofl6 or 7 seconds, this) long pime constant is of little concern. However , 
ey occasional field record has shown a short damped train of 60-Hz 
oscillations, presumably due to power-—demand switching in the relatively 
distant power~line distribution systems, with a resultant step change in 
the amplitude of the local 60-Hz interference field: 

A second limitation of the moditied notch filter 2s Eat eeue 
Noninal22 cipple due to ypulse excitation can be driven to 52 — 102 ae 
the pulse repetition rate is a .precise sub-multiple of 60-Hz. TO -daue’, 
this has not been observed in field recordings, perhaps due to the random 
time distribution of the natural ELF transients and the short coherence 


times of the low-frequency components. 
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APPENDIX B 


THE "GYRATOR" NOTCH FILTER 


B.1 SCHEMATIC AND TRANSFER FUNCTION 


The basic circuit of the so-called "gyrator" notch filter, using 
integrated-circuit operational amplifiers such as the common type ny Gel ees 


issshown in Figure B.l. 


Ry Ro 
Vin | oC 
Ry Amplifier 
WA | a ibs Section 
"Gyrator 
Section 





Figure B.1 The Gy cater Notch Filter 


The voltage transfer function for this notch filter can be 


shown to be: 


+ NG = L R 
go vege ee ae + ws 





[C4C,R,R 
Tey Te. Te ee wales (B-1) 
s2 + se? Gee ae Ww 
VC4CoR,R. 


where Wo, = C{CoR, Re. 
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B.2 POLE-ZERO LOCATIONS 


From the transfer function for H(s), it is evident that a zero 
will occur (the desired frequency notch) if the centre term of the 


numerator is made equal to zero. This will occur if: 


R, + Rg - RoR,/Ry = 


Ry t Rs RgR3/Ry 


die. ce lt 
= Ro, so that the requirement for 


For convenience, make resistors Ry 


a frequency notch becomes simply that: 


Fromenqn.sbolthe filter QO can be written as: 
One yCzCoR Re 
Gj) (Pa + Rye Rs) 
A further simplification results if resistance values are chosen so that 


The transfer function then reduces to: 


W(s) = =p 8 (B-2) 
EWA ory Moby dhryct ares 
where NOW Wo = 1/(RyVC4C2), and the expression for Q becomes: 

; 4 Cy 


This leaves the 0 defined by the ratio of capacitances Cj and Cy, while 


the notch frequency is related to their product via Wo = 1/(Ry{C{C>) - 
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Since the gyrator -section gain is manipulated by the capacitance 
ratio Co/Cy5 the complex-pole location relative to the zero can be 


readily defined in terms of C5/Cy. 






Cy/C1=400 C/C=1000 C/C,=10000 
Q=5 Q=8 0-75 


Figure B.2 Pole-Zero Locations for the Gyrator Notch Filter 


B.3 PRACTICAL LIMITATIONS 


There are two practical considerations which limit the applicat— 


ToOnmmOLetnescyrstom mouchm. 1 Cer: 


(a) The maximum useful Q is about 20. From the expression Repel ap atk 
terms of the capacitance ratio (Eqn. B-3) it can be seen that for Q = 20 
the capacitance ratio C,/Cy = 6,400. With resistor Ry = 100Kohms (about 


the highest practical value in view of the operational amplifier input 


resistance) and f, = 60 Hz, the required capacitance values become: 
C5 = RE uF 


Such a value of Cy can be secured in a metallized-foil capacitor, but it 
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is close to the practical limit of cost and size.(Electrolytic capacitors 
are not acceptable because of the A.C. non-polarized operation, and 
because of their capacitance instability.) Even with a capacitance ratio 
of 6400/1 it is difficult to secure the two units with similar temperat— 
ure coefficients of capacitance. As a result, the notch-frequency drift 


with temperature may be a problem. 


(b) The maximum permitted level of the interference component of V;, 
is a function of Q. This is related to the large capacitance ratio Co/Cy 
which implies that the amplitude of the 60-Hz oscillation (in this case) 
at the output of amplifier A2 (the voltage Vy in Pte Bele must peavey 
much greater than its amplitude at the (+) input of amplifier Al. For a 
maximum signal amplitude of perhaps 10 V (peak) at Vo, the maximum 


amplitude of the 60-Hz component of V;, will be: 


V 10 
Vin(60 Hz) + = 30. = 0.5 V (peak) 


Tt is necessary thatethis interference-level restriction be observed in 
any application of this type of notch filter, since peak clipping of V5 


can introduce transient distortions into the output signal V.. 


B.4 ADVANTAGES 


For applications where a 0 of less than 20 is adequate, the gyrat- 
or notch filter is a simple, stable approach, with relatively few compon- 
ents. Minor frequency adjustments can be easily made if resistor Ry; is 
variable; alternatively, if the required capacitance is sufficiently low, 


capacitor C, can be chosen as a trimmer capacitor. It should be noted 
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that to secure a valid null at the notch frequency, resistor R, may have 
to be somewhat reduced from its theoretical value, to allow for the fin- 
ite input resistances of the operational amplifiers. 

Figure B.3 shows that the pulse response is free of leading and 
trailing-edge differentiation, and the "shock-excited" 60-Hz ripple is 
sufficiently low for most applications. Note that this 60-Hz ripple is 
a constant percentage of the pulse amplitude, and is not a feed through 
of 60-Hz interference from the input signal; in the illustration of 


Fig. B.3(a) no 60-Hz interference was present on the input pulse. 


(a) IBshzeoquare— 
7 Wave Input Signal 


f Notch-Filter 
(b) Output 


Figure B.3 Pulse Response of the Gyrator Noteh Filter 
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